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CALCULATION OF SOUND PROPAGATION IN STRUCTURES 


by L. CREMER, Miinchen 


Summary 


The different types of structure-borne sound waves are described and the wavelengths deduced. 
Especially the radiation into the air and the reciprocal transition are discussed. 
The impedances of simple types of structure are calculated and also the spectrum of impact 


sounds. 


Finally the paper deals with the reduction of longitudinal and flexural wave amplitudes 
afforded by change of section, transition from rod to plate, corners, elastic layers and impeding 


masses. 


Sommaire 


On décrit les différents types d’ondes sonores qui sont transmises par les éléments de cons- 
truction et on en calcule les longueurs d’onde. On étudie spécialement le cas du rayonnement 


dans lair et de la transition en sens inverse. 


On calcule les impédances d’éléments simples de construction ainsi que le spectre des sons 


d’impact. 


Enfin, le présent article traite de la réduction de l’amplitude des ondes longitudinales et de 
flexion qui est produite par les variations de section, le passage d’une tige a une plaque, les angles, 


les couches élastiques et les masses réflectantes. 


Zusammenfassung 


Es werden zunichst die verschiedenen Wellenarten, die bei Baukonstruktionen hauptsichlich 
interessieren, beschrieben und die Wellenlaingen berechnet. Dabei wird besonders die Frage der 
Abstrahlung in Luft und die des reziproken Uberganges aus Luft diskutiert. 

Es werden dann einige Impedanzen berechnet und aus ihnen wird das Kérperschallspektrum 


fiir StoBvorginge abgeleitet. 


SchlieBlich wird die Kérperschalldammung behandelt, die sich bei Longitudinal- und Biege- 
wellen ergibt bei Querschnittswechsel, beim Ubergang Stab— Platte, an Ecken, bei elastischen 
Zwischenlagen und bei aufgesetzten Sperrmassen. 


INTRODUCTION 


A great part of the noise which annoys us in 
buildings or cars is due to the fact that sound 
travels through the supporting or separating cons- 
tructions. This is true for all those noises in par- 
ticular which are directly produced within these 
solids, such as knocking, rubbing, etc. 

In constructional acoustics little attention has 
been paid thus far to the study of the problem 
how the wave propagation actually takes place in 
solids; this may be due to the fact that most of 
the disturbing noises are produced in air and 
finally reach our ear through the air. 

Furthermore, the physics of sound propagation 
in solids is far more complicated than in the case of 
sound propagation in gases and liquids. In the 
case of continua of unlimited extension there are 
two different types of waves, namely the longi- 
tudinal waves (dilatation waves) also occurring in 
liquid or gaseous media and, due to the possibility 
of transmitting shearing stresses, transverse 
shear waves (rotational waves). These two types 
of waves are mostly produced at the same time 


and when one type strikes an obstacle or a 
reflecting surface the other type arises as a 
secondary effect. If the reflecting surfaces lie close 
together, as is the case with all plates and rods, a 
great variety of resultant wave types will be 
produced. 

In general the useful representation of a pro- 
blem of wave propagation in finite structures 
depends on the ratio of their dimensions to the 
respective wavelengths. Since the audio-frequen- 
cies cover 9 octaves, we have to deal with very 
different problems. If low frequencies are of 
interest a construction can be considered as one 
unit, i.e. as a mass or a compliance. In this case 
the whole construction can be treated like a 
coupled system with a finite number of degrees of 
freedom. This treatment does not differ from that 
usual in the mechanics of oscillations. 

On the other hand, if very high frequencies are 
of interest, in particular in the field of ultrasonics, 
all dimensions are large as compared with the 
wavelength and it will be best to use laws of 
propagation proper to a three dimensional con- 
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tinuum. For this case seismology has done a great 
deal of useful preparatory work. 

The characteristic domain of “sound propa- 
gation in constructions” is the medium range. 
Most of the constructional elements, no matter 
whether these are buildings or steel constructions 
of ships, consist either of rods with very different 
cross-sections or of plates. The appropriate repre- 
sentation of problems concerning these structures 
and the frequencies which are of chief interest 
depend upon the fact that, compared with the 
wavelength, they are considered to be small in 
respect to the small lateral dimensions and to be 
large in respect to their longitudinal dimensions. 
Only one or two dimensional propagation prob- 
lems are of interest in so far as sound in con- 
structional elements is concerned. 


I. WAVE-LENGTHS 


1. Longitudinal waves 


If the acoustical properties of a given structure 
are to be determined, first the ratio between the 
dimensions and the wavelengths in question has 
to be considered. If the disturbance occurs in a 
certain narrow frequency range, the problem is 
relatively simple to solve. If the frequency range 
is unknown, the whole audio-frequency range or 
that part of it which is of interest in architectural 
acoustics, i.e. 100 to 3200 c/s, has to be con- 
sidered. 

The computation of wavelengths can only be 
based upon certain physical properties of the 
building materials. This immediately involves a 
first difficulty. The density 9 of most materials is 
a quantity well-known or easy to determine. 
Things are quite different with Young’s modulus 
E, the modulus of rigidity G, and the coefficient 
related to both. The 
results of earlier measurements of Young’s mod- 


of lateral contraction 


ulus mostly refer to its static value and have 
moreover proved to be much dependent on the 
load. 

In architectural acoustics, however, the dy- 
namic value of Young’s modulus has to be con- 
sidered and this may differ considerably from the 
static value. The dynamic value is best deter- 
mined from the longitudinal natural vibrations 
of a free-free rod. Its length / equals half a wave- 
length A, or a multiple of this distance: 

AL 
—on 2° (1) 


From the associated frequency fh 


=n 31 (2) 
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we obtain the velocity of propagation cy of a 
longitudinal wave: 


cr =VE/o. (3) 
Young’s modulus can be calculated from cy if the 
density 0 is known. It is, however, easy to refer 
all the formulae to cz and 0 instead of E. So it is 
sufficient to note the values of cz, which more- 
over represent a more obvious physical con- 
ception. a 
Such determinations of cz, for various building 
materials have been made by R. Scumipr [1]. 
Rods suspended by means of wires were set into 
vibration at one end by a_bone-conduction 
receiver, and the vibrations were received by a 
pick-up at the other end. Every resonance, i.e. 
agreement of driving and natural frequencies, 
resulted in a. narrow peak on the record. Thus he 
was able to determine many values of cy within 
a large frequency range: 


Material cL 

Tron 5170 m/s 
Reinforced concrete 3709 m/s 
Bricks with mortar 2350 m/s 
Timber (fir) lengthwise 4909 m/s 


These values, however, were obtained from single 
specimens of each material. Only the value for 
the metal, i.e. the iron, may also be applicable to 
specimens of other composition. The term “con- 
crete’, even when restricted to heavy concrete, 
involves materials differing in density by a factor 
of 2 and in Young’s modulus by a factor of 2.5. 
As long as we can suppose that both properties of 
different compositions vary in the same sense, the 
value given above would not differ by more than 
12%. For bricks, however, the data given in the 
literature vary considerably, perhaps owing to 
the differences in mortar. E, LUscke [2] e.g. 
gave cy =4300 m/s instead of the above value. 
The data concerning timber are similarly different. 
Compared with the above value obtained by 
Scumipt the well-known pocket-book “Hiitte” 
gives cy = 4180 m/s for fir and 3380 m/s for oak. 
In the direction transverse to the fibres the stiff- 
ness and velocity of propagation may be con- 
siderably lower. For plywood, which is nearest to 
an isotropic material, an average value of 
cy = 3000 m/s can be taken. 

Since the known data are partly contradictory, 
partly not precisely defined, it is to be expected 
that new data gained by systematic tests ac- 
cording to Scumipt’s method will soon be 
published. 

At present we only can refer to the data given 
above when estimating the effects to be expected 
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in architectural structures. In most cases we have 
to restrict ourselves to an estimate of the order 
of magnitude since on the one hand the theoreti- 
cal presuppositions are mostly not fulfilled and 
on the other hand the composition and the dimen- 
sions of the materials and the building methods 
are far from being uniform and exact. 

Under these circumstances some minor differ- 
ences can be disregarded. The difference, for 
example, between the velocity of longitudinal 
waves in a rod, as mentioned above, and that in 
a plate 

er, =VE/[o(1—n*)] (4) 
which with ~ =0.3 is about 5% higher, can be 
ignored in comparision with the uncertainties in 
E and 9. The velocity of longitudinal waves in an 
infinite medium, which is again somewhat higher, 
is of no interest in most buildings. 

The calculation of the wavelength from equa- 
tions (2) and (3) is easy but if it is required re- 
peatedly it is convenient to read it from a dia- 
gram. Fig. 1 shows the wavelengths plotted 
against the frequency as straight lines on a 
logarithmic scale. 
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Fig. 1. Wavelengths of longitudinal waves. 


This diagram shows, for example, that in con- 
crete the wavelength is 37 m at 100 c/s and is 
still 1.15 m at 3200 c/s. From this result it can 
be deduced that a pillar 3 m high can be con- 
sidered as a quasi-stationary element at 100 c/s, 
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i.e. it can be considered as a mass or a spring 
according as its base is soft or rigid. Further, it 
can be deduced that a partition wall h=20 cm thick 
is, even at 3200 c/s, still thin with respect to the 
wavelength so that the necessary restriction for 
the validity of the formula for longitudinal waves 
is not surpassed. This restriction requires about 


Ap> 3h. (5) 


Finally the magnitude of the wavelength is of 
interest for estimating possible resonances. At a 
frequency of about 600 c/s the above mentioned 
pillar would be equal in length to half a wave- 
length, in other words it would vibrate at a natu- 
ral frequency if it is equally terminated at both 
ends either by a hard or a soft medium. 


2. Transverse and torsional waves 


Fig. 2 is a similar diagram showing the wave- 
lengths of transverse waves (shear waves). They 
were derived from the longitudinal waves by put- 
ting “4 = 0.3 for all materials, because more exact 
data were not available. This means: 


Ey ei a (lh -62 es 2a (6) 
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Fig. 2. Wavelengths of transverse waves. 
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These pure transverse waves are of little impor- 
tance in architectural structures. The conditions 
are quite different in vehicles and machines 
especially because at the same time cr is equal to 
the velocity of propagation of torsional waves in 
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structural elements of a circular symmetric cross- 
section. An iron shaft, e.g. 5 m in length, be- 
tween two heavy masses would show its funda- 
mental resonance at 320 c/s. 

In other cross-sections, e.g. in a flat steel strip, 
the torsional waves are propagated more slowly. 
For a narrow rectangle of breadth b and height H 
we have 


er, =) G/g- 2b/H. (7) 
With values of say b=5 and H=50 mm the 


velocity is one fifth of that of the pure transverse 
wave. A bar of this cross-section having a length 
of 1 m would have its fundamental natural fre- 
quency at 320 c/s. 


3. Flexural waves 


Flexural waves are the waves of highest impor- 
tance in architectural acoustics. Their velocity of 
propagation depends not only on the material 
but also on the thickness of the plate or bar and 
the frequency. The dependence on frequency is 
associated with that on thickness since in a sound 
field every length has to be referred to the re- 
spective wavelength, e.g. that of longitudinal 
waves. The ratio h/A, is equivalent to hf/c,. The 
velocity of propagation of flexural waves is given 


by 


cp =e, V1.8hf/er =\1.8cr hf (8) 
and the wavelength of flexural waves by 
dp =\/1.8crh/f. (9) 


The graphs of Ap against f had to be plotted on 
several pages. Fig. 3 refers to iron plates of vary- 
ing thickness. On account of the law of dispersion, 
i.e. the dependence on frequency of cz, the wave- 
lengths in question of one type of plates are 
covered by the ratio Ayog:A3299 = 32:1. Moreover 
we see that the flexural waves are so short that 
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even at low frequencies several wavelengths can 
occur in a structural element. Thus the consider- 
ation of semi-infinite structural elements is jus- 
tified. A sheet of iron 1 mm thick, for example 
has A,o9 =30 cm. A plate of 1.2 x 1.2 m? of this 
material shows at least three wavelengths in 
every direction. 

Fig. 3 is also useful for determining resonances. 
For an iron bar supported at both ends we have 


L=nde/2 - (10) 


which corresponds to equation (1). For a small 
steel plate 5 mm thick and 25 cm long, sup- 
ported at the edges, we obtain a fundamental 
natural frequency of about 190 c/s. If the plate 
is clamped at both ends the natural frequency 
would be 2.25 times higher. The boundary con- 
dition of perfect clamping, however, can hardly 
be fulfilled. For this reason a value between the 
two has to be expected but cannot be determined 
more exactly. 

In the same manner Fig. 4, 5 and 6 refer to 
plates of concrete, brick and plywood and to 
thicknesses h which are of interest. 

The thickness of the structural elements in 
question puts the limit of validity of our repre- 
sentation of the flexural waves within the audible 
range: the validity range is limited by the in- 
equality 


in > 6h. (11) 


Fig. 4 shows this limit and clearly illustrates that 
the present representation is valid for a wall 
20 em thick only up to about 1 ke/s, and for a 
wall 10 cm thick up to nearly 2 ke/s. Similar 
results illustrated in Fig. 5 are obtained with 
brick walls of a thickness equal to half a brick 
(half bat) or a quarter brick (quarter bat). 

The law of dispersion described by equation (8) 


cannot be valid up to arbitrarily high frequencies. 
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Fig. 3. Flexural wavelengths Fig. 4. Flexural wavelengths 
in iron plates. in concrete walls. 


Fig. 5. Flexural wavelengths Fig. 6. Flexural wavelengths 
in brick walls. in plywood panels. 
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This can be concluded from the fact that this 
equation would describe a velocity of propagation 
approaching infinity. The highest velocity of prop- 
agation occurring in solids is, however, that of 
the longitudinal waves (without any lateral con- 
traction). The velocity of the flexural waves could 
at most tend asymptotically towards this value. 

In fact, the velocity of flexural waves tends 
towards a lower limiting value which is (for 
 =0.3) given by 


cr=0.57 cr. (12) 


The subscript R characterizes a special type of 
waves called Rayleigh waves after its discoverer. 
The Rayleigh wave is a surface wave (as illus- 
trated by Fig. 7). This asymptotic behaviour of 
a wall struck by sound waves is understandable. 
At low frequencies, or with a wall thin compared 
with Az, both sides of the wall vibrate in phase. 
Normal to its surface the wall acts as a diaphragm. 
The thicker the wall with respect to A;, i.e. the 
higher the frequency, the smaller is the depth of 
penetration of the forced vibration. 


wee 


Fig. 7. Rayleigh surface waves. 


The velocity of propagation of the Rayleigh 
waves is reached, within about 10%, when 


0:3 Ap. (13) 


With a wall of concrete 40 cm thick this condition 
is attained above about 2.8 ke/s; in other words, 
it is just at the practical limit of thicknesses and 
frequencies occurring in architectural acoustics. 


4. The radiation from a plate into the air 
and the reciprocal transition 


An important part is played by the ratio of the 
lengths of flexural waves to the wavelengths of 
sound in air of the same frequencies. For this 
reason Fig. 3 and the following figures show the 
wavelengths in air by a broken straight line. The 
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slopes of these lines are steeper than that of the 
parallel lines describing the wavelengths of flex- 
ural waves. For this reason there is a critical 
frequency for every thickness of plate, charac- 
terized by the equality of both wavelengths. The 
critical frequency f. can either be seen from the 
points of intersection in Fig. 3 ff. or calculated 
from an equation obtained by putting equal the 
formulae: 


da =V/1-8 cr bf = Ay = col f 
where cy, denotes the velocity of sound in air. 
Hence: 


fe = co[1.8 er h. (14) 


The importance of this critical frequency may 
be understood in regarding the radiation of an 
unlimited flexural wave [3], [6] and [15]. If its 
wavelength Ag is larger than the wavelength in 
air A, corresponding to the same frequency 
HvyYGHENS’ principle proves that there is a wave 
radiated in air on either side of the plate forming 
the angle 

& = arcsin (Ap/Ap) 


between the direction of propagation and the 
perpendicular to the plate. If 7) approaches Az, 
the angle 0 moves towards its maximum value 
of 90°. 

Furthermore if /, surpasses Az the angle ? be- 
comes imaginary which means that radiation fails 
completely. The contrary motions of neighbour- 
ing parts of the plate only produce a quasi-steady 
field of flow in the surrounding medium, shifting 
the particles to and fro without compression. 

With respect to the limited dimensions of the 
partition walls we cannot expect a sudden change 


(15) 


between a frequency range with radiation and 
another without radiation. Yet experiments show 
that the radiation of plates excited at the centre 
increases with increasing frequency and moreover 
they show that this increase is considerable in the 
proximity of the critical frequency. Fig. 8 illus- 
trates the result of such an experiment [4]. A plate 
of aluminium 5 mm in thickness, 1 m in breadth 
and 1.7m in height, was arranged as a partition 
wall between a reverberant and an anechoic cham- 
ber. The plate was excited in the middle at a single 
point by means of the moving coil of a loud- 
speaker. The alternating force could be estimated 
from the intensity of the alternating current and 
thus the supplied power N, could be calculated. 
The power N, radiated from the plate was cal- 
culated from the mean value of sound pressure 
measured at several points in the reverberant 
chamber. Curve a of Fig. 8 shows the quotient 
N,/N, thus obtained (disregarding a constant) 
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Fig. 8. Frequency response of the radiation of flexural 
waves. 


plotted on a logarithmic scale. In spite of the 
relatively small dimensions of the plate we see 
clearly the sudden increase by 10 dB of radiation 
in the proximity of the critical frequency which is 
about 2.6 ke/s. 

In order to gain a second proof of the theory, 
the pressure p, in the anechoic chamber was also 
measured at a distance of 2 m from the driving 
point under an angle of 60° with respect to the 
normal on the plate. The result is represented 
in curve b of Fig. 8 and clearly shows a maxi- 
mum value at the frequency: 

f =fi{smm* 60°. 

According to the law of reciprocity we also may 
expect a better transition from air-borne sound 
in flexural waves of solid plates above the critical 
frequency. Both phenomena are of great impor- 
tance for the by-pass-transition in buildings over 
side-walls and ceilings common to adjacent 
rooms. K.G6sELE! succeeded in checking this 
phenomenon. He used an aluminium plate 2 mm 
in thickness which was placed through a slot of a 
fourfold partition separating two small model 
rooms (Fig. 9a). Both rooms are reverberant so 


‘IT am very much indebted to Herrn K. G6sELE, Stutt- 
gart, for the permission to use his results taken from a 
paper not yet published. The experiments in question were 
subsidised by the ‘“‘Forschungsgemeinschaft Bauen und 
Wohnen, Stuttgart”. 
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that the sound waves in the room containing the 
source struck the plate at all possible angles of 
incidence. In Fig. 9b the solid curve shows the fre- 
quency response of the difference in level between 
both rooms diminished by the tenfold logarithm of 
the absorptive power of room B. We see that above 
4 ke/s this difference suddenly decreases about 
20 dB and rises again above 7 ke/s. The critical 
frequency may be calculated as 6.2 ke/s. 
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Fig. 9. By-pass-transition, measured by GOSELE in model 
chambers. 


(a) Experimental set-up: (1) aluminium plate, 
(2) absorbing material, (3) rubber. 


(b) Solid curve: frequency response of the difference 
in sound level between both rooms, corrected for the 
absorption in room B, 


dotted curve: frequency response of the reduction 
factor of the aluminium plate arranged as a simple 
partition between the two model rooms. 


Furthermore the dotted curve in Fig. 9b re- 
presents the frequency response of the reduction 
factor of the same plate arranged as a simple 
partition between the two model rooms instead 
of the fourfold partition shown in Fig. 9a. In this 
case better transition from air to plate and also 
from plate to air occurs at the same position of 
the plate. 

Regarding formula (14) we may remark that 
a plate radiates the better, the thicker it is. On 
the other hand a thicker plate requires more 
energy to set it in vibration. The above relations 
give an explanation of the fact that impact sound 
caused by hammering a thicker sheet of iron is, in 
spite of lower amplitude, equally as audible as 
sound caused by hammering of thinner sheets 
which show larger amplitudes. 
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Moreover equation (14) explains the fact that 
in water sounds coming from structural elements 
are of higher pitch than sounds radiated by equal 
elements in air. In water the critical frequency is 
higher by the factor (1460/340)? =18, since cy 
denotes the velocity of propagation in the sur- 
rounding medium. Thus in water the radiation 
of the lower partial tones is feebler. 

The relative positions of the wavelength in air 
as against the wavelengths of flexural waves in 
iron plates on the one hand and on the other hand 
in walls and ceilings constructed of bricks or 
concrete are quite different. This is not caused by 
the difference in material but by the quite differ- 
ent order of magnitude of the plate thicknesses. 
The values of cy in masonry are even smaller. In 
almost the whole frequency range the transition 
between air and the wall is rather easy. Only the 
heavy weight and high stiffness counteract an 
observable increase of sound transition as com- 
pared with plates of sheet iron. 

With plywood plates (see Fig. 6) the critical 
frequencies are distributed between about 400 
and 6000 c/s. Therefore, the sound radiating 
properties of thicker wooden plates as well as the 
transition of sound into them can be diminished 
and their sound-insulating properties can be im- 
proved by an artifical reduction of the velocity 
of the flexural waves. This can be accomplished 
by cutting of grooves as illustrated by Fig. 10, 
where the measured reduction factor for a homo- 
geneous plywood plate (1.5 cm) and for a plate 
with grooves is plotted against frequency [5], 
since the velocity of flexural waves in inhomo- 
geneous plates is: 


CB =VBim : V2axf. 
Here m denotes the mass per unit length, or area, 
in the case of a rod or a plate, respectively, and B 
the flexural stiffness, i.e. the quantity given by 
the relation of a static bending moment and the 
resulting curvature, which in the case of a rod is 


given by the product of Young’s modulus and the 
moment of inertia. The critical frequency can be 


Co /% 

rs ie 

By grooving the plates the mass is practically not 

altered, the flexural stiffness, however, is con- 

siderably diminished which results in a higher 
critical frequency. 

This particular measure is an illustration of the 
general rule that it is often advantageous to use 
structural elements of low flexural stiffness if they 
are available. This measure can be applied in a 


(16) 


written: 


(17) 
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building of low sound-insulation in order to 
protect certain rooms from noise propagated in 
the walls or the skeleton of reinforced concrete. 
This aim can be attained to an astonishing degree 
by providing an internal covering of the walls or 
the ceiling by light-weight plates of wood-shav- 
ings which are nailed to laths and coated by 
plaster on the surface directed towards the 
room [7]. The chief reason for this action of the 
additional plates is that they are of low flexural 
stiffness relative to their mass. 
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Fig. 10. Reduction factor of plywood panel without and 
with grooves. 


In the case of inhomogeneous plates or build- 
ing materials of unknown cz, the critical fre- 
quency can be determined in the simplest way by 
a measurement of the sag 7m in the middle of a 
strip of the length | made of the material of the 
plate and supported at its ends. The appropriate 
formula reads: 


fe = 5200 Vm /L2 (18) 


where 7m is measured in centimeters and | in 
metres. This formula says that a strip 2 m in 
length has to show a sag of about 10 cm if the 
critical frequency is to be above 4000 c/s. The 
abovementioned panels or wainscots can be made 
from such plates of low flexural stiffness since 
they are unloaded structural elements. One pre- 
caution only has to be observed, namely, long 
nails must not be driven into the wall at arbitrary 
spots. The application of pictures and other 
decorations, however, will remain unrestricted if 
two horizontal strips 20 cm in height are provided 
for long nails, one just below the ceiling, and the 
other at sight level. All other spots can be reached 
by inserting chains or wires. This demand is likely 
to be fulfilled by every architect, whereas the 
further recommendation to mark the nail-proof 
strips may not be generally adopted. 
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Il. IMPEDANCES 


1. Some simple formulae for impedances 


A second factor of interest in all investigations 
of structure-borne sound is the mechanical imped- 
ance 

7 = Big (19) 
of the respective structural element at the point 
of application of an alternating force F. (v is 
velocity of particles. The bar denotes complex 
quantities.) 

The impedance of an infinite rod set into longi- 
tudinal vibration is 


be = mc, =//m-EA (20) 


i.e. equal to the product of the mass per unit 
length and the velocity of propagation or equal 
to the root of the product of the specific stiffness 
and mass (A area of cross-section). This “input 
impedance” is in the present case identical with 
the mechanical characteristic impedance, i.e. the 
quotient of the force divided by the particle 
velocity in a progressive wave. 

Similarity, for a flexural wave travelling in a rod 
or a plate, the quotient of the lateral force 
divided by the particle velocity is given by the 
product of the mass per unit length and the 
velocity of propagation: 


Zn = mep =m Bo. (21) 


Since for this type of wave, however, the veloc- 
ity increases as the square root of the frequency, 
the impedance of a rod with respect to excitation 
of flexural waves by lateral forces is the higher 
the higher the frequency. 

In a flexural wave the transport of energy is 
represented by the sum of the product of the force 
and the particle velocity and the product of the 
moment M and the angular velocity w, the latter 
being also in phase. It might be suggested that an 
impedance should be introduced differing by the 
square of a length, viz. the impedance with re- 
spect to moments, equal to the quotient of the 
moment divided by the angular velocity in a 
progressive wave. This impedance would read: 

Zou = M/w = Bles =!mB*/o%. —_ (22) 
(A subscript M has been added to point to the 
moment. Correspondingly a subscript Ff should 
be added in equation (21) when confusion is 
possible.) This impedance of moments decreases 
with growing frequency. Thus a rod is the more 
resilient with respect to moments the higher the 
frequency. 
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The geometric mean of the two impedances 
yields a quantity which involves the structural 
data of the rod only: 


Zar: Zem =) Bm. (23) 


It involves only the root of the product of stiff- 
ness and mass as in equation (20). 

If a lateral force acts upon the free end of an 
infinite rod, the quotient of the force and the 
particle velocity is not equal to the expression of 
equation (21). The reason is that a quasi-station- 
ary near-field is also excited. The latter adds to 
the resistive part (real part) of the input imped- 
ance a reactive part (imaginary part) which is 
equal in magnitude and in dependence on fre- 
quency. Moreover, a factor 1/2 turns out to be 


added. Thus we have 


Zz s ha): (24) 
If the point of excitation is in the middle of a 
rod (infinite in both directions) the impedance 


takes a value four times higher: 


Zo = 2mcz(1+i). (25) 


One of the most important and, at the same 
time, unexpected results of the theory of struc- 
ture-borne sound is the fact that the mechanical 
impedance of a plate (infinite in all directions) 
excited at a single point is real and independent of 


frequency [8], [9], [10]. 
Z) =8\) Bm. (26) 
This formula, too, shows the root of the specific 
stiffness and mass, already considered in equation 
(23). Equation (26) can for a homogeneous plate 
of the thickness h be transformed into: 
pines 


— CLoh?. (26a) 

Since the thicknesses of the plates in question 
lie between 0.5 mm and nearly 50 cm, i.e. vary 
by 10°, the values of Z) vary by nearly 10° in spite 
of the small differences in the quantity cp. Fig. 11 
shows the values in question for iron and ply- 
wood, Fig. 12 those for concrete and bricks. The 
thickness h has been plotted as the abscissa. The 
diagram has been subdivided because the values 
for iron and wood lie in a different range to those 
for walls and ceilings. Z, has been plotted as the 
ordinate. The validity of the graphs is based on 
two assumptions. First the wavelength of the 
flexural wave has to exceed six times the value of 
the plate thickness; in other words, the upper 
frequency limits are the same as in the dia- 
grams 3 to 6, Secondly, the point of excitation 
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must be at a sufficient distance from any edge, 
and the attenuation during the propagation and 
at the edges has to be so considerable that any in- 
fluence of reflected waves can be neglected. This 
assumption implies that any focal concentration 
of energy is avoided, in other words, that e.g. a 
circular plate is not excited at its centre. 


2. Spectrum of impact sound 


Knowledge of the mechanical impedance is 
important for several reasons. First, it is the 
basis of an estimation of the force acting upon the 
structural element, or of the particle velocity, 
caused by a transmitter of sound, the internal 
impedance of which can only be estimated 
roughly. Thicker structural elements, the mechan- 
ical impedance of which mostly exceeds the in- 
ternal impedance of the transmitter, require an 
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‘Fig. 11. Mechanical impedance of iron-plates and plywood- 
panels. 


estimation of the force; thinner elements of 
lower impedance that of the particle velocity. 
Secondly, the reaction due to the properties of 
a measuring instrument attached to a structural 
element can be judged. The reaction can only be 
neglected if the mechanical impedance of the in- 
strument is small compared with that shown by 
the object under test at the point of attachment. 
Moreover, the process caused by the impact of 
a mass my striking the structural element at a 
velocity vg, can be calculated, especially its fre- 
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Fig. 12. Mechanical impedance of concrete and brick walls. 


quency spectrum, which is of chief interest. Since 


_at any rate all values of Z in question increase 


with frequency to a lesser degree than does the 
mass impedance of the striking mass (iw my), this 
process can be dealt with as if the striking mass 
were already in contact with the structural ele- 
ment before the impact and experiences a very 
strong but very short pulse containing all fre- 
quencies in the same magnitude. The spectrum 
of the resulting velocity amplitude simply com- 
prises the resulting impedance iwm,+ Z, in the 
denominator, the absolute value of the amplitude 
of a partial tone reading: 


Mo Voo 


[ao] = 


—. (27) 
mt |iwmy+ Zo| 

At high frequencies this expression decreases as 
w—1, At low frequencies and in the case of con- 
stant Z,, i.e. in a rod set into longitudinal vibra- 
tion and a plate set into transverse vibration, we 
obtain a constant value, whereas in a rod set into 
transverse vibration we have a value increasing 
towards lower frequencies as w—™. In the latter 
case the lower frequencies are relatively stronger 
and the impact sound is hollower. Similarly, and 
in accordance with experience, equation (27) 
shows that the decrease as wm! begins later the 
lighter the striking mass m, and the higher Zp, 
i.e. the heavier and stiffer the structural element 
struck by mp». This results in a clearer tone. 
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The noise striking the human ear is moreover 
influenced by the laws of radiation from the 
respective structural elements, and this radiation, 
as already mentioned above, depends on the ratio 
of the length of the flexural waves to the asso- 
ciated wavelengths in air. 

Furthermore, attention has to be called to the 
fact that in measuring spectra by means of octave 
filters, or third-of-an-octave filters, the band 
width of the filters increases with the frequency. 
This means a distortion of the curve to be meas- 
ured because the portion of energy passing 
through any one of the filters is multiplied by f, 
and thus the amplitude of the partial tone com- 
puted by means of equation (27) appears multi- 
plied by f%. The difference, however, of two 
spectra plotted on a logarithmic scale, in other 
words the “improvement” obtained by a certain 
measure, is independent of this distortion. 

The same can be said of the influence on the 
radiation if the measure intended for improve- 
ment does not alter the field of flexural waves 
on the radiating area. 

The simplest measure of improvement consists 
in the insertion of a resilient layer between the 
striking mass m, and the radiating structural 
member which be characterized by Z). Denoting 
the stiffness of this resilient layer by Ky we obtain 
another expression for the resulting impedance 
acted upon by the impact Fi: 

Fy [vo =icom,+ (ze is a 
ve Ly (28) 


This expression, however, cannot be substituted 
for the denominator of the spectrum, for it would 
characterize the spectrum of the particle velocity 
‘vo of the mass my. What we want is the particle 
velocity of the structural element v,, and this is 


related to wo inversely as the associated partial 
impedances: 


Vo] Oe l ie ae 7) | : Lo = ua . (29) 
aut Ze 1+ iw Z,/Ky 


Hence the spectrum in which we are interested 
reads: 


Mp Voo 


‘a = =a Sa 
©" ge \[iomy + Zo/(1 + iv Zy/Ko)] (1 +i Z/Ko)| 
Mp Xo (30) 


x \i@my + Z,(1—w?my/K,)| ' 


Fig. 13 illustrates this expression for constant 
real values of Z,, i.e. specially for a given plate. 
In contrast with the case of direct impact on the 
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Fig. 13. Structure-borne sound spectrum of an impact on 
a plate with resilient layers of stiffness Ky (mp strik- 
ing mass, Vp its velocity, Z) impedance of the plate 
and t = m,/ Zp). 


plate this spectrum decreases as w * at high 
frequencies, and this decrease begins the earlier, 
the lower Ky, i.e. the softer the inserted resilient 
layer. This advantage is compensated by the risk 
of an increase, sometimes even a resonance peak, 
at low frequencies. Since according to our as- 
sumptions losses are ignored, the total impact 
energy is contained in the structure-borne sound 
waves irrespective of the possible insertion of an 
elastic layer. However in the case of an elastic 
layer being used the energy is shifted towards the 
lower frequencies, which is of high importance 
for the audibility of the impact noise, first be- 
cause the energy is shifted towards the frequency 
ranges of feebler radiation, and secondly, chiefly 
because the energy is shifted towards the fre- 
quency ranges of lower audibility’. 

This shifting can be calculated quantitatively 
by means of the present formulae, if the elasticity 
of the intermediate layer is known. Therefore, it 
is recommended to measure the elasticities of all 
impact sound-insulating materials available on 
the market and to have the data published in the 
prospectuses, 

Moreover, there is a certain elasticity in every 
plate of wood, concrete, or even iron, so that the 
idealized spectrum given by equation (27) cannot 
be attained in practice. 


III. REDUCTION FACTORS 


1. Pure longitudinal transition 


The knowledge of the mechanical impedances 
finally permits us in many cases to estimate the 


* A more important measure of improvement is offered 
by the “floating floor” [11]. 
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transition of sound from one structural element 
to another. The transition of longitudinal waves 


is governed by the simple “matching” formula: 


T= (31) 
where r is the generally complex reflection factor, 
i.e. the quotient of the particle velocity of the 
reflected wave and that of the incoming wave, 
Z, the characteristic impedance with respect to 
longitudinal waves in the first medium, e.g. a 
rod, and Z, the input impedance of the second 
medium subsequent to the rod. From r we obtain 
the ratio of the energy which is not reflected but 
transmitted to the second medium divided by the 
incident energy; in other words the so-called 
transmission coefficient: 


d=1—|7? (32) 


which for real values of Z can be transformed 
into 


d = 4/(VZ,/Z_ + VZ,/Z,)°. (32a) 


From d we finally obtain the reduction factor 


RS WOlog 5. (33) 

Considering, for instance, the transition of 
sound from one rod into another of equal material 
but five times bigger in cross-section, we obtain 
from equation (31) to (33): r=—4/6; d=20/36 
and Rw3dB. We see that this considerable 
change in cross-section results in a very small 
insulation effect only. 


2. Transition from rod to plate (sound-bridges) 


Equally simple is the transition of a longitu- 
dinal wave from a rod into a plate normal to the 
rod, although the vibrational motion in the 
second medium consists of a cylindrical flexural 
wave superimposed on a quasi-stationary near- 
field, since the process is governed by Z, only, 
which according to equation (26) has a real value. 
The most important property of the present 
example is the fact that here, too, total trans- 
mission of energy is possible in spite of the dif- 
ferences between the structural elements and 
types of waves occurring. This total transmission 
occurs in the case of perfect matching of the 
mechanical impedances: 


Z,=Z, (34) 
which in the present case means 
m, CL, = —= 05 ¢1, ho. (34a) 


V3 
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In the case of equal media this condition is ful- 
filled if the cross-sectional area A, of the rod is 
made equal to the square of a length about 
1.5 times the thickness of the plate 

4 12 

A, = —=h;z ® (1.5 h,). (35) 

y3 
This dimension is not only easy to design but even 
within the order of magnitude usual in all struc- 
tures involving pillars for support of plates and 
rods between plates. Already from the first 
example we have seen that mismatching even up 
to 1:5 does not result in much insulation. Thus it 
is to be expected that the usual combinations of 
pillars and ceilings in buildings (and in vehicles) 
lead to considerable transmissions of sound and 
noise. It is, therefore, necessary to choose very 
extreme mismatch in order to secure high insula- 
tion. 

This statement is of importance for the studs 
between double partitions, which are already 
known as “sound bridges” on account of this 
disadvantage. These bridges should be so de- 
signed that the characteristic impedance of the 
stud is much smaller than that of the plate or, on 
the contrary, is much higher. In the first case, the 
second plate acts as a rigid terminal for the stud, 
and such a stud usually acts like a spring when it 
is short compared to the longitudinal wavelength. 
In the second case, the stud is not hindered by 
the second plate from vibrating, and thus it acts 
by its inertia. In both cases, calculations show 
that the sound insulation is better, the longer the 
stud, and since all linear extensions in a sound 
field have to be referred to the wavelengths in 
question—in the present case longitudinal wave- 
lengths—the insulation grows with the frequency. 

States of mismatch can be arrived at from two 
directions. Which of them is more suitable de- 
pends on the properties of the plates connected 
by the stud. Between thick walls it would, of 
course, be useless to design connecting studs of 
high mechanical impedance. Such connecting 
studs should be as thin and resilient as permitted 
by their static destination. Even unusual designs 
are worth considering; this is confirmed by 
AsTon’s experiments with his so-called “butter- 
fly-shape ties”, i.e. wires bent into bows for 
higher resiliency [12]. 

In order to connect thin plates, e.g. double 
shells of vehicles, it is mostly easier to choose the 
impedance of the stud to be high in comparison 
with that of the plates, in other words, to connect 
the plates by heavy structural elements. This rule 
seems to be still relatively unknown. Based on 
experience of statics, of thermal conduction, ro 
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of electric current, it has often been concluded 
that the thicker element would result in more 
transmission of sound, but this, however, is err- 
oneous. Already in 1937 E. Meyer [13] had 
shown by experimental tests that the insulation 
of a double partition made of plywood is lower 
with connecting wooden studs than with iron 
connecting elements and gave a qualitative ex- 
planation derived from the better matching be- 
tween the wooden studs and the plywood plates?. 


3. Corners and other junctions 


Equation (31) also permits us to estimate the 
transmission of sound energy entering a rod in the 
form of longitudinal waves and transmitted to 
another rod perpendicular to it. The second rod 
may be connected to the first rod either by a 
swivel link or by a rigid T-joint. The quantity Z, 
has to be taken from equations (24) or (25), 
respectively. In the first case, with rods of equal 
properties, we obtain 


— ocr en(l et), eae 
er dnes( lel) oeaeien st) 


(36) 


where 8 denotes the ratio of the velocities of 
flexural waves and of longitudinal waves. This 
ratio increases as the root of the frequency. The 
new parameter reads 


f = ene = V1.8 hfe 


Generally this quantity is small and according to 
equation (11) limiting the representation of pure 
flexural waves, it is bounded by the inequality 


6 <0.3. (38) 


(37) 


Thus we have a mismatch in any case since the 
input impedance of the second rod set into 
flexural vibrations is always low compared with 
the longitudinal characteristic impedance of the 
first rod. With increasing frequency the matching 
conditions are improved owing to the input im- 
pedance growing as the root of the frequency. 
Equation (36) yields a quantitative expression of 
the transmission coefficient: 


sy 2B 
dpzp =1—|r| ~ 14+ B+ p/2° (39) 


The double subscript LB is to indicate that the 
waves arrive as longitudinal ones and depart as 
flexural waves. 


3 A complete theory of “sound bridges” should also take 
into account that not only forces but also moments can be 
transmitted by the connecting elements. It would be much 
more complicated to deal with this problem exactly. The 
contribution of the moments, however, seems to be in 
general of less importance. 
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According to the reciprocity theorem this 
relation also holds for the energy transmission 
from a flexural wave arriving in the first rod to a 
longitudinal wave departing in the second rod: 


(40) 


The frequency response of dzz described by 
equation (39) is illustrated by a full line in Fig. 14. 

The figure also shows (by a broken line) the 
transmission coefficient dzz calculated for two 
rods or plates connected by a rigid joint. This is 
the case occurring most usually in practice. Un- 
fortunately its calculation is more difficult than 
that of a swivel joint. The six boundary condi- 
tions, viz. the continuity of the two particle 
velocities perpendicular to each other, of the 
angular velocities, of the two forces perpendicular 
to one another and of the moments, are so intri- 
cate that the input impedance Z, is also dependent 
on the data of the first medium and leads to a 
very complicated expression. Nevertheless a rough 
estimate of the question of matching between the 
longitudinal characteristic impedance of the first 
rod— Z;,—and the impedance with respect to 
flexural waves of the second rod— Zpgr2—already 
yields a true picture of the frequency response, 
viz. the transmission described by dzg decreases 
with decreasing frequency. 

In the figure /? has been chosen as the abscissa. 
The respective frequencies depend not only on f 
but also on the material and its thickness. For a 
corner formed by two walls of concrete 10 cm 
thick we have 


dip == dpr . 


ce ak B2 ~ 20000 B? c/s. 


In other words, the highest value permissible in 
this calculation, i.e. 62 =0.1, is associated with 
a frequency of about 2000 c/s. For two plates 
of sheet iron 1 mm thick we obtain 


f =5120002/0.18 c/s. 


0.6 
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Fig. 14. Transmission coefficient versus frequency (para- 
meter /*= 1.8 hf/e,) for a longitudinal wave in a 
rod to a flexural wave in another rod perpendicular 
to the first one. 
solid line: rods connected by a swivel link, 
dotted line: rods connected by a rigid joint. 
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In this example the upper frequency limit of the 
present calculation lies in the ultrasonic range at 


290 ke/s. 


The other limit of the calculated frequency 


response is given by 6%=0.0005. In the first 
example it corresponds to 10 c/s, in the second 
to 1450 c/s. Below this frequency the transition 
of longitudinal waves into flexural waves, and 
vice versa, can in general be neglected, or ex- 
pressed in terms of a boundary condition: each 
of the rods hinders the other from a transverse 
motion. 

In the case of two plates joined at a right angle 
and rigidly connected with respect to flexural 
vibrations an energy transmission by way of the 
product of moment and angular velocity remains 
possible. Thus a flexural wave arriving in the first 
plate can be partially transmitted into the second 
plate. Moreover, it is not surprising that the for- 
mula for dgz involves the ratio of the impedances 
with respect to the moments as defined by 


equation (22): 
2 
Vee VER) 
Cpil/ CB2 CB2}/ CB1 
describing the condition of matching which evi- 
dently is important in the present problem. 


Similarly it is clear that even in the case of per- 
fect matching: 


dzz = 


B,/em = B,/cx: (42) 


i.e. for instance with structural elements of equal 
dimensions and the same material, there is not 
total transmission but we have 


dsp =0.5 (43) 


since a transmission of energy by way of the 
product of force and particle velocity is totally 
hindered by the corner. It is worth mentioning 
that in the case of equal structural elements the 
insulating effect of a corner is only 3 dB. In this 
case, too, the insulation is not much augmented 
by those degrees of mismatch which may usually 
occur. If, for example, the second structural 
element is of double thickness and the same 
material, the ratio of impedances. with respect 
to the moments is augmented by 2’: =5.7:1 but 
the insulation only rises to a value still somewhat 
below 6 dB. 

Similarly worth keeping in mind are the numer- 
ical values describing the transition of energy at 
a T-joint of equal structural elements. The re- 
flected portion is 5/9 instead of 5/10 as above, 
and the transmitted portions are 2/9 each. Finally 
in a joint of four equal structural elements in the 
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form of a cross, the energy of a flexural wave 
arriving in one element is so distributed that 5/8 
are reflected and 1/8 is transmitted in each of the 
three other elements. 

Quantitative confirmations of the present theo- 
retical results are still lacking. We are, however, 
able to state that qualitatively the processes are 
in accordance with the general experience of pro- 
pagation of sound in buildings constructed of one 
material only (e.g. reinforced concrete). No suf- 
ficient insulating effect can be expected from the 
changes in cross-section and direction which are 
mostly given in buildings. 


4. Elastic layers 


An essential improvement in sound insulation* 
can only be obtained by the insertion of special 
structural elements. The insertion of elastic layers 
has proved useful in practice. For longitudinal 
waves the effect is in perfect analogy with the 
insertion of heavy masses between two rooms for 
the insulation of air-borne sound. Formula (44) 
represents the reduction factor valid for the 
present example: 


wh, Ey 


(44) 


where z, is the characteristic impedance and EF, 
Young’s modulus of the interrupted structural 
element (pillar or wall), E, Young’s modulus and 
h, the thickness of the elastic layer, and A, the 
length of longitudinal waves in the medium 1. 
This formula is quite similar to the well-known 
formula of the insulation of air-borne sound by a 
wall of the density 9,: 


One formula shows the influence of the stiffnesses, 
the other that of the densities. In both cases the 
effect is increased by the thickness of the layer for 
any given material, and as the thickness has to 
be referred to the wavelength, it grows with the 
frequency, too. 

If, as an example, we consider a cork layer 
3 cm thick, for which we may assume a modulus 


Ee = 3-107 N/m? (= 3 - 10° dyne/cm?) 


and if this layer is inserted in a pillar or wall of 
concrete, the above approximate formula yields 
an insulation of 7 dB at 100 c/s and 27 dB at 
1000 c/s. The assumptions of equation (44), viz. 
that both parts of the structural element inter- 

4 Recently the term “isolation” has been proposed for 


solid-borne sound, whereas “insulation” is kept for air- 
borne sound. 


330 


rupted by the elastic layer extend over several 
wavelengths, are not always fulfilled. Neverthe- 
less the relatively simple equation (44) is useful 
in permitting us to estimate whether a noticeable 
insulating effect is to be expected at all. 

The insulating effect with respect to flexural 
waves of an elastic layer is also characterized by 
the fact that the reduction factor does not reach 
higher values except at higher frequencies, as is 
illustrated by the example of Fig. 15. The in- 
crease is even steeper and the difference between 
the behaviour of flexural waves as compared with 
longitudinal waves is, moreover, characterized by 
a frequency of total transmission occurring in the 
range of low, but still audible, frequencies. 


10cm cork 


concrete | _ concrete 


3cm 
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Fig. 15. Reduction factor versus frequency for an elastic 
layer of cork between two concrete walls forming 
a straight line. 


1600 3200 6400ck 


This is astonishing since the intermediate layer, 
at any rate, causes a disturbance of the sound 
field, so that an alteration of the field pattern is 
to be expected. This is the case, indeed. It is, 
however, a special feature of the flexural waves 
that such a disturbance consists not only of 
reflected waves but also of a quasi-stationary 
near-field which after being generated does not 
need any further supply of energy. In exceptional 
cases the disturbance in front of the obstructing 
layer consists exclusively of such a near-field. This 
means that no energy is reflected, and thus the 
total incident energy is transmitted’. 


5 Recently this has been verified experimentally by 
M.-L. EXNER and W. Boume [16]. 
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In the present example the associated fre- 
quency of total transmission can be calculated 
approximately from 


Ano = Hh E,|Ep- (46) 


The frequency of total transmission is therefore 
given by: 


——— 2 —— 
Use h oe / E? 
haa gee ee 2° (47) 
wt KE, Hh Fy O01 Hh 


This frequency is lower, the lower the modulus 
and the thicker the layer. It is also dependent on 
the data of the interrupted structural element. 

It is desirable to try to place this frequency at 
a sufficiently low value, not only to avoid the 
“effect of total transmission” but also so that the 
slope of the frequency response, which also is a 
function of the ratio f/f, will begin earlier. The 
response, as illustrated by Fig. 15, can be de- 
scribed by the function: 


R = 10log {1+ 1 —(f/fo)"F) (48) 


which even one octave above the frequency of 
total transmission can already be given approxi- 


mately by 
R30 log (f/f): (48a) 


With flexural waves it is also possible that the 
field behind the obstructing layer consists of a 
quasi-stationary near-field only. In this case the 
obstructing layer is vibrating but does not cause 
any flexural waves in the following medium, i.e. 
it completely hinders the transmission of flexural 
waves. Thus, from the point of view of sound 
insulation, the “‘effect of total insulation” is as 
welcome as its counterpart, the effect of total 
transmission, is undesirable. 

The interruption of a pillar or wall by an inter- 
mediate layer also results in a similar effect of 
total insulation, but this lies in a frequency range 
which is beyond the limits of validity of our cal- 
culation. Nevertheless the approach to this effect 
is already observable in the considerable slope of 
the reduction factor vs. frequency curve accord- 
ing to equation (48a). 

The effect of total insulation as a particular 
feature of the flexural waves can also occur in the 
frequency range in question even at low frequen- 
cies. This is shown by the calculation of the 
equally important problem of an elastic inter- 
mediate layer at a corner, e.g. between a ceiling 
and its supporting walls. This example, dealt with 
and illustrated by Fig. 17, clearly shows the 
physical conditions necessary for the occurrence 
of total insulation. Fig. 16 has been drawn to ex- 
plain this example. The deflections, shown in the 
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figure at intervals of a quarter-period, represent 
a pure near-field and obey the law: 
—2x/A 
160. ss 
This field is caused when the bending moment 
acting at x= 0, 


a) on Ps Ze" 
Miao.B pe 8 (3) sien (49a) 


sinwt. 


and the alternating force acting at the same point, 


oM 2n\* . 
sa aie B (;;) sinwt (49b) 


are in phase (with respect to the signs chosen 
here) and satisfy the proportion 


My: Fy =An:2a. (50) 


n Zig, 


Gi 


Fig. 16. Quasi-steady motion of a rod driven by a bending 
moment and connected with an elastic layer at 
the end. 


In general it is difficult to predict the frequency 
at which these relations are fulfilled. At a corner 
with an elastic layer the relations are exception- 
ally clear since the only quantity given by the 
incident flexural wave is the moment whereas the 
force is generated by the deformation of the 
elastic layer only. Now we have: 


F)=—Kig (51) 


where K is the stiffness of the layer. The force 
resulting from equations (51) and (49) is of equal 
phase to that described by equation (49b) which 
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is required in the case where only a near-field and 
no propagated wave is caused. Both expressions 
are in perfect agreement, i.e. only a near-field is 
caused, when we have: 


K=B a (52) 


AB 


in other words at the frequency of total insulation 
described by 
1 K *Is —i, 


In the arrangement illustrated by Fig. 17 two 
values of K come into question. If the energy is 
transmitted from the x-direction into the y-direc- 
tion, the rigidity 


(52a) 


(53a) 


hinders the vibration of the second structural 
element; for a wave arriving from the y-direction 
the stiffness with respect to pressure 


K,~ E,H|h (53b) 
has to be substituted. By virtue of the reciprocity 
theorem the insulating effect is the same in either 


direction. For this reason we obtain two frequen- 
cies of total insulation, viz. 


oe 
E\-o1 
AOS c7 Ve) aK (54a) 
according to equation (53a) and 
pea a aie eed 
1a! 


according to equation (53b). The example illus- 
trated by Fig. 17 is based on the same data as 
Fig. 15, viz. a 3 cm cork layer between concrete 
plates 10 cm thick. A comparison of the formula 
(54) with the formally similar expression in 
equation (47) shows that the two frequencies of 
total insulation lie somewhat below and some- 
what above the frequency of total transmission 
shown in Fig. 15. Hence both frequencies lie in 
the lower frequency range. 

With a perfect total insulation R ought to tend 
towards infinity at the above frequencies. The 
height of the insulation peaks, however, is re- 
stricted by the generation of secondary longitu- 
dinal waves and by the internal losses in the cork 
layer. 

Assuming that in cork 10 % of the deformation 
energy is irreversibly dissipated, we may expect 
insulation values of 50 dB. 

In the range of high frequencies the insulation 
is also satisfactory since here the slope of the 
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insulating effect of an elastic layer with growing 
frequency is favourable. The frequency response 
of the reduction factor tends to the values given 
by the approximate formula: 


(55) 


with cork 


BS 


0 50 100 200 400 800 1600 cls 


Fig. 17. Reduction factor versus frequency for a corner 


formed by concrete walls. 


full line: with 
dotted line: without 


an elastic layer of cork. 
This expression is similar to that given in equa- 
tion (44) for the insulation of longitudinal waves 
by an elastic layer, only differing in that equation 
(44) involved the wavelength of longitudinal 
waves, which is inversely proportional to the fre- 
quency. In the present example we have the 
wavelength of the flexural wave which is inver- 
sely proportional to the root of the frequency, and 
this results in a smaller increase of R with fre- 
quency according to 


R = 10log f + 10 log (5.5 h? H— Ej’ 0)" Ex?). (55a) 


In the example of Fig. 17 the second term of this 
sum proves to be so large that in the interesting 
range the associated straight line always runs 
above 36 dB. 

For the present problem of flexural waves a 
simple solution has been found that is similar to 
the corresponding problem of longitudinal waves. 
The reason is that at high frequencies the trans- 
mission due to the product of force and particle 
velocity is no longer effective since the impedance 
due to the compliance of the elastic layer de- 
creases with w' and is finally so small in com- 
parison with the increasing impedance Zpr; that 
the plates can be regarded as being free with 
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respect to transverse vibrations. The difference 
between the impedance with respect to the mo- 
ments of the flexural wave Zp which decreases 
with m—” and the corresponding impedance of an 
elastic layer of the thickness h is not so high. For 
a plate of the thickness H the former impedance is 


By fem: ~ E, H3/12 on 
and the absolute value of the latter is 
|M, |(w,—w,))= E, H3/12ho. 


Moreover, it is obvious that the reduction factor 
depends on the ratio of these two quantities just 
as in equation (44) it depends on the ratio of the 
longitudinal characteristic impedance 2,4 and 
the longitudinal impedance AE,/wh, due to the 
compliance. The fact that we have only to sub- 
stitute one of these ratios for the other without 
any alteration of the numerical factor is rather 
difficult to derive but easy to keep in mind. 

In the above calculations the elastic layer has 
been dealt with as pure compliance. This involves 
the assumption that the thickness of the layer is 
always small compared with the lengths of longi- 
tudinal waves in the material in question. The 
velocity of longitudinal waves in cork is about 
the same as in air. Hence a 3 cm layer can be re- 
garded as small compared with the wavelengths 
of the frequencies of interest in architectural 
acoustics. This condition is not fulfilled with 
sponge rubber, in which the appropriate velocity 
can be reduced to 60 m/s. Moreover, it is usual to 
apply thicker layers of this material. Hence it is 
possible that already at medium frequencies such 
a layer must be dealt with as a continuous me- 
dium. The thickness of a 5 em layer of a material 
of this low velocity would at 600 c/s equal half 
a wavelength and permit total transmission if the 
damping could be ignored. But owing to internal 
damping, we need not be afraid of total trans- 
mission. Nevertheless the fact that the insulation 
can be considerably lower than would result from 
the relations and formulae given above has to be 
considered. 

The usefulness of applying intermediate layers 
of cork between ceilings and supporting walls has 
been well recognized for a long time, but only in 
a qualitative sense. Their quantitative caleul- 
ation would be an advance in the art of building. 
The present formulae are intended to aid in solv- 
ing this problem. Moreover the effect of total 
insulation derived in this paper could be used in 
principle. 

Considering the effect of mismatch of the im- 
pedance of the intermediate layer and that of the 
interrupted medium, we tend to insert layers of 
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high acoustical “softness”. In other words, the 
insertion of layers of little stiffness and big thick- 
ness between solid bodies promises a successful 


insulation. Where no static forces have to be | 


transmitted an intermediate layer of air is the 
best material. Nevertheless such a layer has a 
certain stiffness which is often undervalued. At 
high frequencies the air particles cannot give way 
by moving sideways. 

For this reason the dynamical stiffness of the 
elastic layer below a “‘floating”’ floor is nearly 
equal to that of the air volume if this is filled 
by a very yielding and porous material such as 


glass wool [4], [11], [14]. 


5. Impeding masses 


On the other hand it is in principle also pos- 
sible, just as in the corresponding problem of air- 
borne sound, to obtain a certain insulation by the 
insertion of materials of higher acoustic imped- 
ance. These materials act as inert masses, and 
their effect can be calculated for longitudinal 
waves from equation (45) taken from the cor- 
responding problem in air, or more exactly from 
the original equation: 

2 
R= 10log [i+ eee 


1 


(56) 


From this formula we can deduce that for example 
an iron plate 3 cm thick inserted in a vertical 
wall of concrete would even at 10000 c/s result 
in a reduction factor of 3 dB only. This small 
effect shows that acoustically there is so little 
difference between iron and concrete that it is 
permissible to deal with reinforced concrete as a 
homogeneous medium. ; 

The insulating effect of the inserted mass can 
be augmented by giving it a larger cross-section 
than that of the interrupted structural element. 
In house building this measure would require 
very heavy and bulky structures for obtaining 
useful insulating effects, whereas in steel struc- 
tures, e.g. in ships, this possibility may be taken 
into account. 

In shipbuilding any insertion of elastic layers 
seems to be almost impossible for fear of leakage. 

Experience has shown that ribs provided to 
give stiffness to the hull and acting as heavy 
masses from the view point of structural acoustics 
result in a high sound-insulating effect. This 
could hardly be expected from equation (56). 
Here, too, the ‘“‘effect of total insulation” of the 
flexural waves causes the high reduction factors 
with respect to these waves. It was this problem 
which led to the theoretical discovery of this 
effect, confirmed by experiments. 
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In this problem the frequencies of total insu- 
lation can be estimated approximately by means 
of the very simple rule that the wavelength of the 
blocked flexural wave is equal to the circum- 
ference of a circle whose radius equals the radius 
of gyration of the additional heavy mass: 


[hip == PA GbR 


(57) 


This simple formula is more difficult to derive 
than equations (54). Nevertheless it is clear that 
equation (50), being the relation’ between the 
moment and the force, involves the ratio of the 
impedance of the hindering object with respect 
to moments (wmyr?) divided by the translatory 
inertia wm, and this ratio is equal to the square of 
the radius of gyration r (my denoting the additional 
mass provided at a rod or, in the case of a plate, that 
of a transverse rib expressed by its mass per unit 
width). At any rate this consideration shows that 
the additional mass must be compelled to follow 
the rotational vibrations of the plate or rod as a 
rigid body. Therefore the radius of gyration can- 
not be obtained by designing a lumped mass 
maintained at a certain distance by spokes. An 
absolutely rigid design is required. The theory as 
well as tests have shown that the length of the 
connection between the additional mass and the 
extended structural element (in the direction of 
propagation) has to be small compared with the 
wavelength. 

A design combining high stiffness with a large 
radius of gyration and a small connecting length 
is that of a wedge fixed at its edge to the struc- 
tural element in which the flexural waves are 
propagated. Its radius of gyration is 0.7 of the 
“height” of the wedge. Its height depends on the 
material, the thickness of the carrying structural 
element, and the desired frequency of total insu- 
lation. For plates of the thickness h equation (57) 
yields: 


(eee ee 
r= 5-\/1.8erh/f (58) 
In the special case of iron we have 
h{em| 
r ~ 150 ee cm]. 58a 
ley] case 


As an example a sheet iron 0.5 cm thick and a 
frequency f= 200 c/s lead to a radius r=7.5 cm, 
and a “height” of the wedge of 11 cm. These 
dimensions are easily applicable in practice. The 
choice of 200 c/s as the frequency of total insu- 
lation involves a useful insulation over the whole 
frequency range in question, the frequency of 
total transmission which occurs in this example, 
too, is mostly two octaves lower so that one 
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octave lower the insulation may be expected to be 
still fairly good. 

At frequencies above that of total insulation 
the insulation value decreases but soon turns to 
a new increase in an asymptotic approach to the 
approximate response we are going to derive. 
Also in this example the range of higher fre- 
quencies is governed by a formula similar to 
equation (45) obtained when considering the at- 
tenuation of longitudinal waves by a heavy mass, 
VizZ.: 


(59) 


R= 20log (F = 

ABm 

where m refers to the mass of the structural 
element carrying the flexural waves, i.e. in the 
case of a rod m denotes the mass per unit length, 
and in the case of a plate the mass per unit area. 
Similarly, as in equation (55), the transmission of 
energy only depends on one of the products. In 
the present case the product of the moment and 
the angular velocity is ineffective, since the im- 
pedance with respect to rotation of the addi- 
tional mass due to its inertia grows as the fre- 
quency w whereas the impedance with respect to 
moments of the flexural wave Zgm decreases with 
w—’?: in other words the last mentioned imped- 
ance is so small in comparison with the first 
mentioned that any angular velocity vanishes. 
Instead of this equation (59), there is involved 
the ratio of the impedance according to trans- 
latory inertia wm, divided by the impedance with 
respect to forces of the flexural waves Zar= mcs. 
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With cp growing as the root of the frequency, 
R also in this example increases with 10logf: 


R = 10log f+ 10log (1.6m m_* B-"*). (59a) 
In Fig. 18 the full line shows the frequency 


response of the reduction factor due to an addi- 
tional mass. The typical features can be clearly 
seen: a frequency of total transmission in the 
range of low frequencies, a peak at the frequency 
of total insulation, and finally the transition to 
an asymptotic slope. The curve refers to the 
structure indicated by the sectional drawing, an 
additional mass of 115 g attached to a flat steel 
rod 2.7 mm x 2 cm in cross-section. The re- 
duction factor of this arrangement could be de- 
termined experimentally. The measured values 
have been plotted and prove that the effect of 
total insulation can be practically brought about 
and its frequency chosen by a prior calculation. 


6. Problems of oblique incidences 


All formulae discussed as yet have concerned 
either rods or plates at normal incidence of the 
waves on the interface.. Since in practice also 
oblique incidence is of importance, its most im- 
portant cases have to be dealt with. The process 
of oblique incidence is governed by a general rule 
known from the physics of optical refraction. 
This rule deals with the transition from a medium 
1 with a velocity of propagation c, into a medium 
2 with a velocity c, and says that the angles of 
incident and outgoing direction are related by 


eee ol: 


Feat 60 
sind, sind, (09) 
This may be interpreted as the equality of the 
“trace velocities”. If the velocity in the second 
medium exceeds that in the first, c,>c,, the 
angle 3, reaches 90° at an angle of incidence 


Bed k- 
d, = arcsin (2) 


2 (61) 
and at higher values of #, no waves pass into the 
second medium. As an example we consider two 
plates of sheet iron, one being of double thickness 
compared with the other. The velocities of flex- 
ural waves are as |/2:1. At an angle of incidence 
exceeding 45° no flexural wave enters the second 
sheet iron. 

But also in the case of equal thickness of the 
two sheet iron plates, i.e. in the case of matching, 
the transmission coefficient depends on the angle 
of incidence. It decreases with oblique incidence 
according to 


(62) 


d= , cos? i. 
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A discussion of further examples which mostly 
involve very complicated problems will be re- 
nounced, since this essay is meant to show that 
many fundamental problems of structure-borne 
sound can be judged and estimated by means of 
simple considerations and formulae. 
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UBER DEN GEHORSINN VON STECHMUCKEN 


Von H. TISCHNER 


Physikalisches Institut der Universitat Tiibingen 


Zusammenfassung 


Messungen des Schallpotentials an miannlichen Tieren einer Stechmiickenart, Anopheles 
subpictus, ergaben, da die Tiere in einem schmalen Frequenzbereich von 150---550 Hz hoéren. 
Die gréBte Empfindlichkeit wird bei 380 Hz erreicht. Sie entspricht der Empfindlichkeit des 
menschlichen Ohres bei 2 kHz. Der Hérumfang vom Schwellwert bis zur Sattigung betrigt 
28 Phon. Die beobachtete BandpaBwirkung ist auf die Abstimmung der GeiBeln zuriickzufiihren. 
Das Fluggeréusch der Weibchen enthalt einen subjektiv deutlich wahrnehmbaren Anteil von 
380 Hz, so dafi dieser Ton aus dem Geréusch der Umgebung herausgehoben wird. Das eigene 
Fluggeriusch wird, da es héher liegt, von den GeiBeln des Mannchens nur wenig erfabt. Die 
beobachtete frequenzverdoppelnde Wirkung des Johnston-Organs geht aus seinem Aufbau her- 
vor. Die Grundfrequenz-der auf die Stifte wirkenden Spangenbewegung ist bei einer GeiBel- 
bewegung in Richtung der Antennenachse starker vertreten als bei einer Bewegung senkrecht 
zu ihr. Das Tier ist trotz der gegeniiber der Wellenlinge sehr kleinen Ausdehnung der GeiBeln 
in der Lage, die Richtung zu bestimmen, aus der das weibliche Fluggeraéusch kommt. Eine be- 
sondere biologische Bedeutung kommt der Uberlagerung zweier Tone zu, die als. Schwebung in 
der unmittelbaren Umgebung von 380 Hz ablaufen. Die relativ grof8e Schwebungslaufzeit des 
Schallpotentials hat physikalische Ursachen, die durch die Resonanzeigenschaften der Antenne 


bedingt sind. 


Summary 


Measurements of the sound potential in male animals of a species of mosquitoes, Anopheles 
subpictus, show that these insects are hearing in a small range of frequency between 150 and 
550 c/s. The maximum sensitivity occurs at 380 c/s. It corresponds to the sensitivity of the 
human ear at 2 ke/s. The difference between the threshold level and the level of saturation is 
28 phon. The noise produced by the flying female mosquitoes contains a clearly perceptible tone 
of 380 c/s so that this tone is of great use to the male insect in spite of the general noise 
of the environment. The noise of flying of the male mosquito will practically not be picked up 
by its own feelers as the loudest tone of the spectrum is at the upper limit of its auditory range. 
The observed effect of doubling the frequency of the Johnston Organ may be understood from 
its construction. The first harmonic of the movement of the clasps working on the pins is larger 
at a motion of the feeler in the direction of the axis of the antennae than at a motion perpendicular 
to it. In spite of the small dimensions of the feeler compared with the wavelength the mosquito 
is able to find the direction from which the female noise of flying is coming. A particular biological 
importance belongs to the superposition of two sinusoidal sounds which form beats near 380 c/s. 
The relatively high time of delay of these beats observed in the sound potential has physical 


origins in the resonance effect of the antennae. 
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Sommaire 

Des mesures faites sur le sens de l’ouie de sujets males d’une certaine espéce de moustiques, 
Anopheles subpictus, ont montré que ces insectes entendent dans une gamme de fréquences 
étroite, de 150 & 550 Hz, la fréquence de plus grande sensibilité étant 380 Hz. Cette sensibilité 
correspond 4 celle de l’oreille de ! homme pour 2 kHz. Il y a 28 phones entre le seuil d’audition et 
la saturation. L’effet de bande passante observé est di a l’accord des flagella. Le bruit du vol de 
la femelle contient une composante nettement perceptible subjectivement a 380 Hz, de sorte 
que ce son ressort du bruit ambiant. Le bruit de vol lui-méme, de fréquence plus élevée, n’est 
presque pas percu par le male. L’effet observé de doublement de la fréquence par lorgane de 
Johnston est dé a la constitution propre de cet organe. La fréquence fondamentale du mouvement 
de serrage agissant sur les tiges est plus marquée quand le mouvement des flagella a lieu dans la 
direction de l’axe que dans la direction perpendiculaire. Malgré la trés-faible dimension des 
flagella devant la longueur d’onde, l’insecte peut déterminer la direction d’ot vient le bruit du 
vol de la femelle. La superposition de deux sons qui forment un battement au voisinage de 380 Hz 
a une importance biologique particuliére. La réponse relativement retardée du sens de louie 
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a ce battement a une origine physique liée aux propriétés de résonance des antennes. 


1. Einleitung 


Die Frage, ob Stechmiicken (Moskitos) héren, 
bildet seit langer Zeit den Gegenstand der Er- 
érterungen von Entomologen. Seit dem Jahre 
1855, in dem C. Jounston [1] das nach ihm be- 
nannte Gehérorgan von Anopheles beschrieb und 
C. M. Cuitp [2] in seiner Leipziger Dissertation 
dieses Organ im Jahre 1894 schon recht genau 
dargestellt hatte, wurde diese Frage gelegentlich 
immer wieder behandelt. Man beobachtete, dai 
Schwarme von Stechmiicken durch Violin-, Cello- 
oder Zitherspiel [3] und durch Singen angelockt 
wurden und als die elektrotechnischen Erzeug- 
nisse EKingang ins praktische Leben fanden, stellte 
man fest, das der von bestimmten Telegraphen- 
apparaten und einem ténenden Lichtbogen abge- 
strahlte Schall Moskitos anzog. Sir H. Maxim [4] 
bemerkte beim Bau einer Kraftanlage, daB eine 
schadhafte Lampe, welche als Telephon wirkte, 
von einem Mickenschwarm umgeben war. 

Untersuchungen von A.M. Mayer [5] mit 
Stimmgabeln aus dem Jahre 1874 zeigten, daB 
die GeiBelenden von Moskitos besonders kraftig 
bei 520 Hz ausschlugen. Bei zwélf Tieren war 
die Resonanzfrequenz ungefahr gleich. Eine klare 
Vorstellung von der Wirkungsweise des Johnston- 
Organs fehlte jedoch bisher. Die fiihrenden ame- 
rikanischen Entomologen L. C. Howarp, H. G. 
Dyer und F, Knap [6] faBten ihre Meinung im 
Jahre 1912 in dem folgenden Satz zusammen: 
Wenn der Schall mit den Antennen empfangen 
wird, wie es auf Grund des Baues der Antennen 
der Fall zu sein scheint, so ist noch zu bezweifeln, 
ob das Tier wirklich auf Schall reagiert oder ob 
die GeiBeln nicht vielmehr auf den Luftstrom 
ansprechen. Es ist auch unklar, wie die Antennen 
auf Schall verschiedener Tonhéhe reagieren. 

Die ersten zuverlassigen Untersuchungen iiber 
das Verhalten von méannlichen Tieren einer 
Moskitoart, Aédes aegypti, des Gelbfieberiiber- 
tragers, wurden von L. M. Rorn [7] im Jahre 1948 


mitgeteilt. Er arbeitete ebenfalls mit Stimm- 
gabeln. Der Frequenzbereich, auf den Aédes 
aegypti bei sinusférmigem Schall anspricht, er- 
weitert sich etwas mit dem Alter der Tiere. Er 
liegt nach den auf statistischem Wege gefundenen 
Ergebnissen zwischen 200---500 Hz. AnschlieBend 
registrierten M. C. Kann und W. OFFENHAUSER 
jun. [8] das Fluggerausch weiblicher Aédes 
aegypti auf einer Schallplatte. Sie fanden bei der 
Wiedergabe mit einem Lautsprecher, das mann- 
liche Insekten der gleichen Art beinahe sofort 
nach dem Einschalten auf den Lautsprecher zu- 
flogen. 

Fir die Beurteilung des Hérvermégens hat 
man von drei verschiedenen Seiten auszugehen: 
Einmal ist die Morphologie des Gehérorgans zu 
studieren, das bei mannlichen Culiciden ungleich 
besser entwickelt ist als bei den Weibchen, zum 
anderen ist das Verhalten der Tiere gegeniiber 
dem weiblichen Fluggerausch zu untersuchen und 
schlieBlich sind mit elektrophysiologischen Me- 
thoden die elektrischen Potentiale zu erfassen, 
die entstehen, wenn Stechmiicken sich im Schall- 
feld aufhalten. Die morphologische Seite wurde 
jiingst von H. Risier [9], ausgehend von der 
Arbeit von C. M. Curip [2], wieder aufgegriffen, 
der im AnschluB an die hier wiedergegebenen 
Beobachtungen die Anatomie des Gehérorgans 
von Anopheles stephensi studierte. Ergebnisse 
iiber das Verhalten der Tiere liegen von L. M. 
Rot sowie M. C. Kann und W. OrFENHAUSER 
jun. vor. 

In der vorliegenden Arbeit wurde der dritte 
Weg beschritten. Mit den Untersuchungen von 
E.G. Wever und C. W. Bray [10] iiber die 
Cochlearpotentiale an der Schnecke von Wirbel- 
tieren als Grundlage wird iiber Messungen der 
Schallspannung an einer Moskitoart, Anopheles 
subpictus, berichtet, die wahrend eines Aufent- 
haltes in Indien am Indian Institute of Techno- 
logy zur Verfiigung stand und die mit A. ste- 
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phensi eng verwandt ist. Die drei Wege miissen 
zusammenfiihren und eine einheitliche Vorstel- 
lung von der Wirkungsweise des Gehérorgans von 
Stechmiicken ergeben. 


2. Elektrophysiologische Untersuchungen 


In Abb. 1 ist die Schaltung zur Messung des 
Schallpotentials an der lebenden Miicke darge- 
stellt. Zwei spitze Stahlelektroden fiihren vom 
Tier zu dem Eingang eines erdunsymmetrischen, 
einstufigen Verstarkers, der sich mit dem Moskito, 
den Elektroden und deren Haltevorrichtung in 


\ 
Schwebungs - 
summer 
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Abb. 1. Schaltung zur Messung des Schallpotentials. 


einem gerdumigen, doppelt geschirmten Fara- 
dayschen Kafig befindet. An den Ausgang des 
zweiten, einstufigen Verstarkers sind handels- 
ubliche Réhrenvoltmeter R.V. und zwei Katho- 
denstrahloszillographen [11] K.O. I und K.O. II 
angeschlossen. Vor dem Faradayschen KAfig ist 
ein Lautsprecher L angeordnet, in dessen Schall- 
feld das Tier liegt. Zur Messung des Schalldruckes 
befindet sich in seiner unmittelbaren Nahe (Ab- 
stand des Moskitos kleiner als ein Sechstel der 
Wellenlange) ein dynamisches Mikrophon M. 
Zwischen 100---800 Hz war die gemessene Span- 
nung praktisch unabhangig von der Frequenz. 
Die an den Eingangsklemmen des ersten Ver- 
starkers auftretende Rauschspannung betrug 
6 wV. Da sich dieser Verstarker im Schallfeld des 
Lautsprechers befindet, und eine Klingfrequenz 
bei 750 Hz auftrat, wurde er in zwei gut schall- 
dicht voneinander getrennte Kasten eingebaut. 

Die Tiere wurden mit Chloroformdampf 
schwach betaubt, mit ihren Beinen auf einen 
Objekttriger geklebt und eine Elektrode in das 
linke Auge gestochen. Die andere Elektrode 
wurde langsam in die Gelenkfalte an der Pedi- 
cellusbasis der rechten Antenne gefiihrt und in 
diese eingestochen, bis das Braunsche Rohr eine 
Spannung anzeigte. Auch Versuche, bei denen 
eine Elektrode in das linke Auge, die andere in 
den Hinterleib des Tieres gestochen wurden, ver- 
liefen erfolgreich, wenn auch die gemessene Span- 
nung viel kleiner war als in der unmittelbaren 
Umgebung des Johnston-Organs. Belegte man 
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die zweite Elektrode mit einer feuchten Graphit- 
schicht, so wurde ebenfalls eine Spannung be- 
obachtet, wenn die Elektrode zart auf den Pedi- 
cellus gelegt wurde. Sie hob sich gerade vom Stér- 
pegel ab. Das Kriterium der ,,richtigen“* Mef- 
stelle ist so ein rein elektrisches und kein opti- 
sches. Es ist méglich, daf der umfangreiche 
Nervenstrang vom Johnston-Organ zum Gehirn 
des Tieres nicht getroffen wurde, da die nach- 
folgend beschriebenen Ergebnisse nicht das zei- 
gen, was von physiologischer Seite als Nerven- 
aktionspotential (Impulsfrequenzmodulation) be- 
zeichnet wird. 

Ersetzte man den Kopf der Micke durch einen 
Leimtropfen der gleichen GréBe und stach beide 
Elektroden in ihn hinein, so blieb der Stérpegel 
auch im relativ starken Schallfeld unabhangig 
von der GréfBe des Schallfeldes konstant. So war 
es sicher, dafi die gemessenen Erscheinungen 
wirklich zum Tier gehéren. 


100uV. 


= aaa f-20% | 


 ——— 
a. 


Storpegel des 
Verstarkers(6uV) 


oa! 


600 Stonpegel des 
fonpegel des 
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P 


i) 200 400 


Abb. 2. Verlauf des Schallpotentials in Abhangigkeit vom 
Schalldruck. 


Das Ergebnis einer MefBreihe zeigt Abb. 2. 
Mit steigendem Schalldruck P, in willkirlichen 
Einheiten (mV) gemessen, erhebt sich das Schall- 
potential wu tiber den durch den Anschluf des 
Tieres erhéhten Stérpegel von 16 2V zuniachst 
linear, um bei einem bestimmten Schalldruck 
einen Sattigungswert zu erreichen. Da kein ge- 
eichter Schalldruckmesser zur Verfiigung stand, 
wurde bei verschiedenen sinusférmigen Ténen der 
Schwellwert des Schalldruckes in der Nahe des 
Moskitos subjektiv bestimmt und tiber den in 
Spannungseinheiten geeichten Schwebungssum- 
mer und die bekannte Fletscher-Munson-Kurve 
auf die absolute GréBe des Schalldruckes, bei dem 
die Sattigung erreicht ist, geschlossen. Dieser 
Sattigungsdruck (Schnittpunkt der gestrichelten 
Linien) hangt stark von der Frequenz ab. In 
Abb. 3 ist er als Funktion der Frequenz darge- 
stellt. Wahrend im Bereich um 150-:-600 Hz der 
Sattigungsdruck relativ grofs ist, zeigt er bei 
380 Hz ein scharf ausgepragtes Minimum. 

Ein Versuch, den Schwellwert zu ermitteln, 
fihrte zu der zweiten Kurve fiir A. subpictus. 


338 
2Harmonische | Grundton 
1008 | | 
| 
~ 80 7) 
g, Sathgung 
Qa 
S 60 
8 
§ 40 
: ~ i 30 Phon 
| 20 menschl. Ohr 
Schwellwert 0 Phon 
oe 1 I a 
40 6080100 200 400 600 1000 2000 4000 100Q0H2 
if 2 
Abb. 3. Empfindlichkeitskurve des Gehérorgans von A. 


subpictus. 


Das Schallpotential wurde oszillographisch kon- 
trolliert und derjenige Schalldruck bestimmt, bei 
dem sich dieses Potential gerade vom Stérpegel 
abhob. Bei 240 Hz gelang dies noch bei dem glei- 
chen Tier, mit dem die Sattigungskurve erhalten 
wurde. Die Messungen erfolgten im Februar 1952 
und mufBten aus 4uBeren Griinden fiir acht Tage 
unterbrochen werden. Nach der Wiederaufnahme 
der Versuche war das Leben im dicht benach- 
barten Dschungel so laut geworden, da Versuche, 
den Schwellwert erneut auf der alten Grundlage 
zu bestimmen, erfolglos verliefen. Der Stérpegel, 
vorwiegend von Insekten hervorgerufen, war auf 
den dreifachen Betrag angestiegen. Ein schall- 
toter Raum stand nicht zur Verfiigung. Der Ver- 
lauf der Schwellkurve wird sich jedoch nicht 
wesentlich von dem der Sattigungskurve unter- 
scheiden. Es ergibt sich so ein Hérumfang von 
28 Phon. Der Schwellwert ist im Bereich der 
gréBten Empfindlichkeit sehr klein. Er ist etwa 
gleich dem Schwellwert des menschlichen Ohres 
bei 2 kHz. Das Mannchen der Stechmiicke hért 
damit das relativ leise Fluggerausch des Weib- 
chens besonders gut, aus dem ein Ton von 380 Hz 
stark hervortritt. 

Uber die Héhe des subjektiv wahrgenommenen 
Flugtones von Culiciden liegen Untersuchungen 
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Abb. 4. Subjektiv wahrgenommene Flugténe ven Stech- 
miicken. 
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von QO. Soravatta [12] (Abb. 4) vor. Die Ton- 
héhe des Fliigelschlages ist von der Temperatur 
abhangig. Wichtig ist, daB bei allen Stechmiicken- 
arten der Unterschied zwischen dem minnlichen 
und weiblichen Flugton deutlich in Erscheinung 
tritt. Die bei 28°C durchgefiihrten eigenen Mes- 
sungen zeigen dies ebenfalls. Der Unterschied 
betragt auch bei A. subpictus etwa 200 Hz. Die 
Flugtonhéhe des Mannchens ist ungefahr 540 Hz; 
sie liegt also an der Grenze seiner eigenen Wahr- 
nehmungsfahigkeit. So ist die Méglichkeit gege- 
ben, daB die mannlichen Tiere, die in den Abend- 
und Morgenstunden in Schwarmen [13] auftreten, 
ihr eigenes Fluggerausch von dem vorbeifliegender 
oder in den Schwarm eindringender Weibchen 
trennen kénnen. Aus den Beobachtungen von 
A.M. Mayer [5] schlieBen wir, daB die hohe 
Empfindlichkeit bei 380 Hz durch die mecha- 
nische Resonanz der beiden GeiBeln bedingt ist. 

Ein Blick auf den zeitlichen Verlauf des Schall- 
potentials liefert Aufschliisse iiber den Umwand- 
lungsprozeB des Schallvorganges im Johnston- 
Organ. Der elektrische Stérpegel bei einer sehr 
kleinen Amplitude des duBeren Stérgerausches 
gibt eine unregelmaBige Kurve (Abb. 5a). Dieses 
und die nachfolgenden Oszillogramme zeigen als 
untere Kurve jeweils den Verlauf des Schall- 
druckes oder eine Zeitmarke. Ahnlich ist das Er- 
gebnis bei einem Schallvorgang mit 140 Hz 
(Abb. 5b). Die Aufnahmen wurden mit einem 
Zweistrahloszillographen und einer Leica- Kamera 
gemacht. Oszillograph und Kamera befanden sich 
in etwa 1 m Abstand von dem Mikrophon. Durch 
das Offnen des Schlitzverschlusses der Kamera 
entstand ein Gerausch, dessen Spektrum zum 
Teil in den Hérbarkeitsbereich der Micke fallt. 
Die Auswirkung dieses Gerausches auf das Schall- 
potential ist in dem rechten Teil des Oszillo- 
gramms zu erkennen. 

Bei einer Besprechung des Tieres mit 240 Hz 
(Abb. 5c) erscheint ein véllig regelmaBiger Ver- 
lauf des Schallpotentials, bei dem die zweite 
Harmonische stark hervortritt. Dies ist charakte- 
ristisch fiir jeden sinusférmigen Vorgang unter- 
halb der Frequenz gréSter Empfindlichkeit. Es 
kommt auch vor, daB bei dicht benachbarten 
Frequenzen der Phasenwinkel des Grundtones 
sich schnell gegeniiber dem der zweiten Harmo- 
nischen dreht. Oberhalb 380 Hz tritt der Grund- 
ton immer mehr hervor; bei 400 Hz ist er fast 
ausschlieBlich vertreten (Abb. 5d). 

Es ist lehrreich, die Reaktion von Stechmiicken 
auf die menschliche Stimme zu studieren. Auf 
Singen des Vokals ,,a“* durch eine mannliche 
Stimme reagiert das Tier mit einer Grundfrequenz 
von 320 Hz und guter Frequenzverdopplung 
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Abb. 5. (a) Stérpegel des Schallpotentials. 
(b) Schallpotential bei 140 Hz. 
(c) Schallpotential bei 240 Hz. 
(d) Schallpotential bei 400 Hz. 
(e) Schallpotential des Vokals ,,a‘‘. 
(£) Schallpotential des Lautes ,,sch“‘. 
(g) Schallpotential des Einschaltstosses eines Laut- 
sprechers. 
(h) Schallpotential des Dschungelgeriusches. 
(i) Schwebung des Schallpotentials 340/391 Hz. 
(k) Schwebung des Schallpotentials 370/392 Hz. 


(Abb. 5e). S-Laute werden wegen ihrer hohen 
Formantbereiche zwischen 5-:-10 kHz nicht mehr 
wahrgenommen; dagegen wird der Laut ,,sch“ 
(Abb. 5f) sehr gut vom Tier registriert. Im 
Johnston-Organ entsteht eine Frequenz von 
770=2 x 385 Hz. Man kann sich das Vergniigen 
machen, die an der Stechmiicke abgegriffene 
Spannung iiber Verstarker einem Kopfhérer zu- 
zufiihren und die Sprache abzuhéren, wie sie im 
Johnston-Organ erscheint. Sie ist nur noch sehr 
unvollkommen verstindlich. Das dem engen 
DurchlaBbereich entsprechende Frequenzband 
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tritt deutlich hervor. Die Anordnung wirkt als 
ausgesprochener BandpaB. Unterhalb der opti- 
malen DurchlaBfrequenz wird also die Frequenz 
verdoppelt, oberhalb bleibt das Potential sinus- 
formig. 

Die BandpaBbwirkung kommt auch bei der Auf- 
nahme von EinschaltstéBen des Lautsprechers 
zum Ausdruck (Abb. 5g). Aus dem unregelmafi- 
gen Verlauf des Schallvorganges wird die Frequenz 
720= 2 x 360 Hz heraussortiert. Wichtig ist die 
Reaktion des Moskitos auf das Dschungelge- 
rausch. Die Messung erfolgte abends 10 Uhr, als 
die Insektenwelt in voller Bewegung war. Auch 
aus diesem Gerausch (Abb. 5h) wird die Frequenz 
2 x 360 Hz ausgewahlt. Die untere Kurve ist wie 
bei den Abb. 5e und f eine Zeitmarke. 

Die Natur achtet streng auf die Trennung der 
Arten. Da die Flugtonhéhe nach den Beobach- 
tungen von QO. Soravatra [12] bei allen Stech- 
mickenarten praktisch die gleiche ist und auch 
zeitlich konstant bleibt, miissen. andere Unter- 
scheidungsmerkmale fiir das Fluggerausch der 
einzelnen Arten vorliegen. Es liegt nahe, die Un- 
terschiede in dem Modulationsrhythmus und der 
Hiillkurve des von den Fliigeln abgestrahlten 
Schallvorganges zu suchen. W. LoTTERMOSER [14] 
beobachtete die grobe RegelmaBigkeit der von 
Waldgrillen und Heuschrecken hervorgebrachten 
Locklaute. Die einzelnen Arten unterscheiden 
sich durch die Form der Hillkurve und die Tast- 
frequenz. Auch zwischen A. subpictus und A. ste- 
phensi war der subjektiv gewonnene Eindruck 
des Fluggerausches bei etwa gleicher Flugton- 
héhe verschieden. Neben der Amplitudenmodu- 
lation ist auch eine Frequenzmodulation denkbar. 

Um einen Anhaltspunkt fiir diese Auffassung 
zu gewinnen, wurden durch die Uberlagerung 
zweier sinusférmiger Spannungen ein sinusférmig 
amplitudenmodulierter Schallvorgang hergestellt 
und das zugehérige Schallpotential gemessen. So 
wurden zwei Spannungen von 340 und 391 Hz 
zwei Schwebungssummern entnommen und ad- 
diert (Abb. 5i). Auch hier tritt im Schallpotential 
die doppelte Frequenz der resultierenden Schwin- 
gung von 2x 365 Hz klar hervor. Die Schwe- 
bungen weisen eine relativ groBe Verzégerungszeit 
von 10,95 ms auf. Véllig anders ist das Bild, wenn 
zwei Spannungen von 370 und 392 Hz iiberlagert 
werden (Abb. 5k). An die Stelle der Schwebungen 
des Schallpotentials treten scharfe Spitzen. Diese 
Erscheinung laBt darauf schlieBen, dab jetzt eine 
besondere Art der Schallumwandlung vor sich 
geht, die fiir den Anlockungsmechanismus von 
Bedeutung ist. Die Einstellung der beiden Fre- 
quenzen ist ziemlich kritisch. Die auch hier be- 
obachtete grobe Schwebungslaufzeit von 18,4 ms 
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ist nicht auf die endliche Fortpflanzungsgeschwin- 
digkeit der Schallpotentiale in den Nerven zu- 
rickzufiihren; sie hat, wie im Anhang dargelegt 
ist, rein physikalische Ursachen, die durch die 
Resonanzeigenschaften der Geifeln mit ihrer Kin- 
spannung bedingt sind. 


3. Zur Deutung der Ergebnisse 


Aus den elektrophysiologischen Beobachtungen 
ergab sich die Empfindlichkeitskurve und die 
frequenzverdoppelnde Wirkung des Johnston- 
Organs. An Hand der von C. M. Cutzp [2] und 
H. Rister [9] durchgefiithrten morphologischen 
Studien lassen sich nun Vorstellungen von der 
Wirkungsweise des Gehérorgans von Stech- 
miicken entwickeln. Diese vor allem beim mann- 
lichen Tier relativ groBen Organe befinden sich in 
den am Kopf vor den Augen gelegenen Pedicellen 


der Fithler (Abb. 6). In ihnen wurzeln die beim 


Abb. 6. Dorsalsansicht des Kopfes von A. subpictus. 


Mannchen besonders lang und dicht behaarten 
GeiBeln. Ein Langsschnitt durch den Pedicellus 
(Abb. 7) zeigt den nicht ganz einfachen Aufbau 
des Johnston-Organs. In der aus dieser Abbil- 
dung abgeleiteten halbschematischen Darstellung 
(Abb. 8) [15] wurzelt die GeiBel G, die im Innern 
ein BlutgefaB und Nerven enth4lt, in einer Fuf- 
platte Pl, die an Spangen Sp verankert ist. Die 
durch die Bewegung der Geibel an der Platte ent- 
stehenden Druckschwankungen wirken iiber die 
ringférmig angeordneten Spangen auf die End- 
fasern F der Stifte St, im auBeren Ring, welche 
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Abb. 7. Langsschnitt durch die Antennenbasis (sagittal). 
Phasenkontrastaufnahme. 


ventral «— dorsal 


Abb. 8. Rekonstruktion von Scapus, Pedicellus und GeiBel. 
Rechts ein Teil der inneren Pedicelluswand und die 
Stifte des inneren Ringes entfernt. Schnittflachen 
wei. 
G Geibel, Pl FuBplatte, Sp Spangen, F Endfasern, 
St, Stifte des iuBeren Ringes, St, Stifte des inneren 
Ringes, Sz Sinneszellen, P Pedicellus, Se Seapus. 
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in Gestalt einer Hohlkugel um die Platte herum 
zusammengefaBt sind. Die ihnen zugeordneten 
Sinneszellen Sz liegen nahe der AuBenwand des 
Pedicellus P. Der ringférmige Scapus Se und der 
Pedicellus werden durch Muskeln bewegt, um der 
gesamten Antenne eine bestimmte Richtung rela- 
tiv zum Kopf zu geben. Neben den duferen 
Stiften St, besteht ein zweiter innerer Ring von 
Stiften St,, die mit ihren Endfasern an der Basis 


der Spangen haften. Ein Querschnitt (Abb. 9) 


Abb. 9. Querschnitt durch den Pedicellus in Héhe der FufB- 
platte. 


Sp Spangen, F Endfasern, Sz Sinneszellen, Pl Fub- 
platte, St, Stifte des iuBeren Ringes. 


durch den Pedicellus in Héhe der FuGplatte Pl 
1aBt die Spangen Sp mit den Endfasern F, die 
Stifte St, und die Sinneszellen Sz erkennen. Das 
Organ ist im wesentlichen rotationssymmetrisch 
um die Antennenachse herum aufgebaut. Eine 


Abb. 10. Vereinfachte 
Abb. 7 und 8. 


Darstellung des 


Lingsschnittes 


H. TISCHNER: UBER DEN GEHORSINN VON STECHMUCKEN 


341 


weitere Schematisierung zeigt Abb. 10. Die Stifte 
St, miinden in die Sinneszellen Sz, die an dem von 
der Platte Pl entfernten Ende angeordnet sind. 

Im Schallfeld verursacht eine Schnellekompo- 
nente in Richtung der Geibel eine Bewegung der 
Platte parallel zu sich selbst. Beim Niedergang 
driicken die Spangen Sp auf die Stifte St, der 
unteren Sinneszellen Sz, bei der Aufwartsbewe- 
gung auf die Stifte St, der oberen Zellen. Ein 
sinusférmiger Ton hat so ein Schallpotential der 
doppelten Frequenz zur Folge, wie es unterhalb 
380 Hz beobachtet wurde (Abb. 5c). Dariiber 
hinaus erscheint der Grundton, der durch die 
Bewegung der Stifte St, hervorgebracht wird. 
Wegen der geringeren Zahl der Stifte ist auch 
seine Amplitude kleiner. In Abb. 5c betragt die 
erste Harmonische 60% der zweiten. 

Tritt eine Schnellenkomponente senkrecht zur 
Antennenrichtung auf, so wird die Platte gekippt. 
Auch die Stifte St, geben nunmehr die doppelte 
Frequenz des Schallvorganges. Daraus folgt, dah 
das Tier die Richtung, aus der der Schall kommt, 
feststellen kann. Der Grundton verschwindet um 
so mehr, je kleiner die Schnellenkomponente in 
Richtung der GeiBel ist. Die Wellenlange des 
Schalles bei 380 Hz ist etwa 1m, so da®B die 
Phasenverschiebung zwischen den beiden Geibel- 
bewegungen keine Rolle spielt. Der Abstand 
zwischen den beiden Geifeln an ihrer Einmiin- 
dung in die Johnston-Organe betragt 0,2 mm. 
Diese Phasenverschiebung ist nur bei gréBeren 
Abstanden der beiden Gehérorgane (Wirbeltiere) 
von wesentlicher Bedeutung. Die Natur hat mit 
der Konstruktion des Johnston-Organs die Még- 
lichkeit geschaffen, sehr kleine Empfanger durch 
die Erfassung des Klirrfaktors richtungsempfind- 
lich zu machen. Dies gelingt um so besser, als 
durch die BandpaBwirkung der GeiBeln aus dem 
weiblichen Fluggerausch nur ein praktisch sinus- 
formiger Anteil von ungefahr 380 Hz fiir die 
Signaliibermittlung abgesondert wird. 

Bei héheren Frequenzen als 380 Hz tritt die 
zweite Harmonische zuriick (Abb. 5d). Es scheint, 
als ob das Johnston-Organ nicht mehr in der 
Lage ist, den raschen Schwankungen des Schall- 
feldes zu folgen. 


ANHANG 


Die zeitliche Verzogerung der Bewegung 
in Resonanzschaltungen 
bei Schwebungsvorgangen 


Die GeiBeln der Stechmiicken sind resonanz- 
fahige Gebilde. Die Masse wird durch die Geibel 
mit ihren Harchen und die Platte (Abb. 8) dar- 
gestellt; die Einspannung der Fufplatte ist der 
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elastische Teil der schwingenden Anordnung. Die 
Bewegung x der Geibel gehorcht, wenn man an- 
nimmt, da sie als Ganzes schwingt, der Differen- 
tialgleichung 


(1a) 
wobei m die Masse, r die Dampfungskraft je 
Schnelleeinheit, ¢ die Riickstellkraft je Weg- 
einheit, P,) die Amplitude der erregenden Kraft, 
w ihre Kreisfrequenz bei sinusférmigem Antrieb 
und ¢ die laufende Zeitkoordinate bedeuten. Die 
Lésung der Gleichung fiir den eingeschwungenen 
Zustand ist 


tl a 
mx+rx+cx = Pye”, 


Aol i(@t—e) 
X= XE 


(1b) 
mit x) als Bewegungsamplitude und ¢ als Nach- 
eilwinkel der Bewegung gegeniiber der Kraft. 
Durch Einsetzen der Lésung in Gl. (la) erhalt 
man mit den Abkiirzungen n=o)/m/c; D= 
r/2Vme im Anschlu8 an K. Kiorrer [16]: 


Le Pe; 
c(1—7?+i2.Dn) ° 


(le) 


Kae 
Mit cx,)/P,)=y ensteht 
1 
; e=arctg ua st kd) 
Jn + 4D? 17 
Die reduzierte kinetische EinfluBzahl y stellt als 
Funktion der normierten Frequenz 7 die bekannte 


Resonanzkurve (Abb. 11) dar; der Nacheilwinkel 
e steigt von Null iiber 7/2 bei 7=1 auf z. 


Af = 


Abb. 11. Verlauf der kinetischen EinfluBzahl y und des 
Nacheilwinkels ¢ als Funktion der normierten 
Frequenz 7. 


Wirken zwei Krafte mit den Amplituden P,, 
und P,, mit den Kreisfrequenzen @, und w, auf 
das schwingende System, so iiberlagern sich 
wegen des durch die Gl. (la) ausgedriickten line- 
aren Verhaltens die EinfluBzahlen y, und y, mit 
den entsprechenden Nacheilwinkeln ¢, und ¢, zu 
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ye “=y,e “+y,e ° (2a) 
Um zu einer einfachen Diskussion der Ergebnisse 
zu gelangen, wird angesetzt, da Py), = Po, ist und 
daB die Krafte den Phasenwinkel Null zuein- 
ander im Zeitpunkt t= 0 haben. So wird die an- 
treibende Kraft in reeller Schreibweise: 


Ff (©) = Po (sina,t+ sinw,t). (2b) 


Wir wahlen die Vorgange aus, deren Einflub- 
zahlen y, und y, gleich-sind, legen uns also auf 
bestimmte Werte 7, und 7, bzw. w, und @, fest. 
Damit wird in reeller Schreibweise der zeitliche 


Ablauf der Bewegung: 
&(t)=y1 [sin (wt —e,) + sin (w2t—e,)], 
oder 


(@1 — 2) t + (e2—&) ‘ 


g(t) =2y, c0s 5 


(2c) 


(1 + Wy) t — (eg + &) 
: 


+ sin 


Abb. 12. Die Schwebungslaufzeit. 


Die Funktionen f(t) und g(t) stellen amplituden- 
modulierte Schwingungen (Abb. 12) dar. Die 
Nullstellen der Hiillkurven liegen bei 


(@,— 9) Ty 


yt = (2n—1)5 (n=1,2,3...) (8a) 


bzw. 
(3b) 


(@,;—@y) T,+ (€—&) =(2n 5 (n=1,2,3 walla 


2 


So ist eine zeitliche Verzégerung 
& 1 &—& 


Wo No-Ny 


Tas Tact (Ae ee 


Wy—W, 


(3e) 


der Bewegung gegeniiber der antreibenden Kraft 
zu verzeichnen, wobei wo7,;=@, bzw. Wotjz= 
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gesetzt ist. Da stets e,—¢,<a ist, wird die 
Schwebungslaufzeit 


1 
INSP Se ae 
2 fo(42—m) 


Sie ist um so gréBer, je kleiner der Unterschied 
zwischen den erregenden Frequenzen ist. Fiir den 
Fall der Abb. 5i ist f>=380 Hz; 7,—7,= 0,135 
und AT=10,95 ms. GemaB Gl. (3d) ist ATS 
9,65 ms zu erwarten. Bei dem Versuch nach 
Abb. 5k wird unter der Annahme, da der Impuls 
dem Maximum der Schwebung entspricht, mit 
fo= 380 Hz; 7,.—7,=0,058 und AT=18,4 ms. 
Die maximal mégliche Verzégerungszeit betragt 
22,6 ms. Die gemessenen, auf den ersten Blick 
iiberraschend groBen Verzégerungszeiten haben 
damit eine eindeutige physikalische Ursache, die 
sich zwanglos aus dem Verhalten des Empfangs- 
apparates als resonanzfahiges System ergibt. Die 
Tatsache, das die gemessenen Laufzeiten nahe 
an den maximal méglichen liegen, deutet auf eine 
scharf ausgepragte Resonanz mit einem steilen 
Anstieg des Nacheilwinkels hin. 

Dem Direktor des Indian Institute of Techno- 
logy, Kharagpur (Indien), Herrn Dr. J. C. Guosu, 
danke ich fiir das Interesse, das er dieser Arbeit 
entgegenbrachte, ebenso den Herren M. Kunpu 
M. Se. und P. C. Duar, die mich bei den Messun- 
gen unterstiitzten. Mein besonderer Dank gilt den 
Herren Prof. Hans H. WEBER (Physiologisches 
Institut der Universitat Tiibingen) und Prof. 
HERMANN WEBER (Zoologisches Institut der 
Universitat Tubingen), die mir in Unterhaltungen 
Gelegenheit zu einer Urteilsbildung auf einem 


(3d) 


‘Aeronautical Acoustics, in particular Jet Noise” 


With this title a symposium was held, under the joint 
auspices of the British Acoustics Group and the Royal 
Aeronautical Society, at the Institution of Mechanical 
Engineers on May 2Ist, 1953. The first paper by Prof. 
E. J. Ricwarps discussed the work being done in several 
Universities on sounds of supersonic jets. He was followed 
by Mr. J. B. GrEATREX who described the measurements 
by Rolls-Royce, Ltd., on the distribution of noise round jet 
engines and its reduction. Mr. G. M. Litiey discussed the 
theory of aerodynamic noise, in particular that which re- 
sults from the interaction of turbulence with shock waves 
and from shear in boundary layers. 


H. TISCHNER: UBER DEN GEHORSINN VON STECHMUCKEN 
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Gebiet gaben, das mir bei Beginn der Arbeiten 


neu war. 
(Eingegangen am 20. Marz 1953.) 
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The remaining papers were concerned with the noise 
associated with actual aircraft in flight; Mr. WARREN’s on 
the ‘bangs’ heard by listeners on the ground as the air- 
craft breaks through the ‘sound barrier’, Mr. HAyHURST’s 
on the sounds heard in the neighbourhood of the starting- 
up and Mr. FLEMING on the spectra and loudness of the 
steady noise produced at the ground by aircraft passing 
overhead. 

A long and interesting general discussion followed rang- 
ing from the theory of the generation of supersonic noise 
to what constitutes a noise nuisance and the legal aspects 
of the problems of its mitigation. 

The papers presented in the symposium will be publis- 
hed by the Royal Aeronautical Society, 4 Hamilton Place, 
London, W. 1. E.G.R. 


THE RELATION BETWEEN THE SONE AND 
PHON SCALES OF LOUDNESS 


by D. W. ROBINSON 
National Physical Laboratory, Teddington, Middlesex, England 


Summary 


Several previous investigators have studied the relation between loudness level and loudness 
sensation but have obtained discordant results. It is thought that a cause of these differences is 
that results are influenced by the test procedure. A new set of experiments is described in which 
a variety of sounds having known loudness level-intensity relations was used to make independent 
estimates of the phon-sone scale law over a wide intensity range. The uniqueness of this relation, 
which is prerequisite to its existence, is strongly demonstrated by the agreement of results for 
entirely different sounds and also by the consistency of the data for widely different test ratios. 
An unexplained difference remains between the scales for sounds ordered in ascending and des- 
cending loudness respectively, as has been found previously. The experiments were conducted in 
a manner directly related to the practical use of the sone scale, and were made with statistical 
control. The nature of the sone scale near the absolute threshold is discussed. A review of previous 
work leads to a single set of data representing tolerably well the results of various authors which 
may serve as a basis for future discussion of the standardisation of the phon-sone relation. 


Sommaire 


Plusieurs chercheurs ont déja étudié la relation existant entre le niveau sonore subjectif en 
phones et l’impression sonore subjective en sones, mais ne sont pas arrivés a des résultats con- 
cordants. L’auteur estime que ce désaccord est di a la différence des méthodes de mesure. Il 
décrit une nouvelle série d’expériences effectuées sur tout une variété de sons, dont on connaissait 
la relation entre l’intensité physique et l’intensité subjective, afin d’établir indépendamment la 
relation entre l’échelle des phones et celle des sones sur une gamme étendue d’intensités. L’unicité 
de cette relation, condition nécessaire de son existence, est nettement démontrée par le bon accord 
entre des résultats relatifs 4 des sons entiérement différents et aussi par la consistance des résultats 
concernant des rapports trés différents. I] reste une différence inexpliquée entre les échelles a 
intensités croissantes et décroissantes, respectivement, différence déja constatée antérieurement. 
On a fait les essais en se plagant dans des conditions yoisines des conditions d’emploi de échelle 
des sones, et sous contréle statistique. On discute la nature de l’échelle des sones prés du seuil 
absolu. L’examen des travaux antérieurs conduit 4 une série unique de valeurs qui représente 
suffisamment bien ensemble des résultats des différents auteurs et qui peut servir de base pour 
une discussion ultérieure en vue de la normalisation de la relation entre les phones et les sones. 


Zusammenfassung 


Bei mehreren friiheren Untersuchungen der Beziehung zwischen Lautstirke und Lautheits- 
empfindung ergaben sich voneinander abweichende Resultate. Es wird angenommen, dal} eine 
Ursache hierfiir in einer Beeinflussung der Ergebnisse durch den Testvorgang liegen kénnte. In 
einer hier beschriebenen neuen Versuchsreihe wurde eine Anzahl verschiedener Schallvorginge 
mit bekannter Lautstarke-Intensitiats-Beziehung benutzt, um in einem grofen Intensitatsbereich 
unabhingige Bestimmungen des Gesetzes der Phon-Sone-Skala durchzufiihren. Die Allgemein- 
giiltigkeit dieser Beziehung, die eine Vorbedingung fiir deren Existenz ist, wird durch die Uber- 
einstimmung der Ergebnisse fiir véllig verschiedene Schallvorginge und auch durch die Konstanz 
der Werte bei sehr verschiedenen Testverhiltnissen tiberzeugend nachgewiesen. Es bleibt eine 
unerklarte Abweichung zwischen den Werten, wenn man die Schallvorginge nach zu- oder ab- 
nehmender Lautstirke anordnet, was schon friher festgestellt wurde. Die Versuche wurden mit 
statistischer Kontrolle und in direkter Beziehung zur praktischen Benutzung der Sone-Skala 
durchgefiihrt, deren Verlauf in der Nahe der Hérschwelle diskutiert wird. Ein Uberblick iiber 
friihere Untersuchungen fiihrt zu einem einzigen Satz von Werten, der ziemlich gut die Ergebnisse 
verschiedener Autoren wiedergibt und als Grundlage fiir weitere Diskussionen iiber die Normierung 
der Phon-Sone-Beziehung dienen kann. 


1. Introduction sound pressure level of an equally loud pure tone 


In recent years it has become common practice of frequency 1000 c/s [1]. The tentative inter- 
to express the loudness of noise by means of the national agreement on the phon scale reached in 
“phon” scale, that is to say, in terms of the 1937 has brought about a welcome uniformity in 
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the terminology of noise measurement. However, 
loudness values expressed on this scale do not 
immediately suggest the actual magnitude of the 
sensation of loudness and this sometimes leads 
to difficulty in interpretation on the part of en- 
gineers and other users who need to translate the 
phon values of noise into terms of common 
experience. For instance, the sequence 10—40— 
70—100 phons suggests a quantity rising by 
equal increments, but these values do not in fact 
correspond to sensations having this simple rela- 
tionship. 

It has often been proposed that to meet this 
difficulty there should be a supplementary scale, 
of which the “sone” scale has been proposed as 
an example [2], by means of which a number 
representing the average observer’s judgment of 
sensation is assigned to a noise as an aid to the 
interpretation of the loudness level. If, for ins- 
tance, the value of one sound on the sone scale 
were double that of another, the first sound would 
be judged by the average listener to be about 
“twice as loud” as the second, and similarly for 
other ratios. 

It was the purpose of the present investigation 
to provide data which may assist in establishing 
the scale relation between phons and sones and 
to provide a basis for any future discussion on the 
standardisation of the scale for practical pur- 
poses. 


2. General 


A number of reports of similar experiments 
have already appeared in the literature [3]---[8]. 
Rather large discrepancies are found between the 
results obtained by different workers, and it 
seems that details of procedure of the subjective 
measurements involved may have a relatively 
large influence on the judgments of observers. In 
this new series of subjective loudness experiments 
an effort has been made to avoid as far as possible 
the causes which are thought to influence the 
judgment of subjects, and in addition the experi- 
mental plan has been so devised as to reveal its 
own self-consistency or otherwise in a rather 
sensitive manner. 

At least three methods of determining the 
relationship between intensity and loudness have 
been described, namely: 


(1) by direct subjective estimation of the inten- 
sity increment required to produce an assigned 
ratio of loudnesses, 


(2) by observing the intensity increment required 
to maintain constant loudness on changing 
from binaural to monaural listening, 
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(3) by comparing the intensity of a complex 
sound built up of a number of equally loud 
components with that of one of the component 
tones judged to be of the same loudness as the 
complex sound. 


Whereas the first method is direct, the last two 
methods, an account of which has been given by 
FLETCHER and Munson [9] are indirect, and their 
validity rests on assumptions about the mecha- 
nism of hearing. There seems to be some justifi- 
cation for considering experiments of types (2) 
and (3) as being more of the nature of tests of the 
assumptions than as a basis for the subjective 
loudness scale itself. 

It was concluded by GrIcEeR and FIRESTONE [5] 
on the basis of direct determinations by method. 
(1) that the concordance of the three methods was 
established. However, in order to reach this con- 
clusion it was necessary to find fault with the 
determinations of previous workers [3], [4]. Since 
the paper by GercEer and FIRESTONE, other deter- 
minations have been reported, notably by 
CuurcHEerR, Kine and Davies [6], Potiack [7] 
and GARNER [8]. Some of these later results show 
more affinity with those of the workers criticised 
than with GrIcER and FIRESTONE’s own results, 
and their conclusion on the equivalence of the 
direct and indirect methods does not now appear 
so certain. 

It therefore seemed desirable to attempt a new 
determination of the intensity-loudness relation 
by a direct method under conditions which were 
as far as possible free from the objections which 
have been or might be levelled against previous 
experimental techniques. 

In practice, the loudness scale has its main 
application in assessing the effectiveness of noise 
reduction projects and other cases involving the 
immediate comparison of two noises according 
to a numbered scale. This suggests that the me- 
thod by which the subjective measurements ought 
to be made is by direct comparison of pairs of 
stimuli, the loudness ratio of which is judged in 
relation to an assigned numeral. This direct me- 
thod also permits proper comparison to be made 
with the results of earlier investigations similarly 
carried out. 

With the exception of Ham and ParkINson’s 
work and a few supplementary measurements by 
Larrp, each of the previous direct loudness-scale 
determinations has been made using earphones. 
There are some limitations attached to this pro- 
cedure. The field of practical importance of noise 
measurement does not ordinarily include sounds 
heard through earphones; consequently it is pre- 
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ferable in laboratory tests to work with free 
sound fields so that in this respect at least the 
conditions of the laboratory tests may approxi- 
mate to those met with in practice. Further, the 
use of free sound fields admits of more accurate 
specification and simpler measurement of the test 
sounds, especially if these are complex in cha- 
racter. 

With regard to the nature of the subjective 
judgments and the possibility of their expression 
in numerical ratios, there exist divergent opinions. 
The standpoint adopted in this paper is that the 
issue may properly be decided by experimental 
investigations. The self-consistency of the experi- 
mental results may be expected to indicate, given 
sufficiently crucial tests, the extent to which the 

resulting loudness scale is real and useful. It 
appears from the present study that experimental 
evidence strongly indicates that the loudness 
scale is real, even though it may still be imper- 
fectly known, and that it can provide a useful 
numerical guide in the consideration of noise 
problems. 

In spite of the scepticism sometimes heard, 
neither the author, nor apparently most other 
experimenters, have found that subjects have 
much difficulty in making repeatable direct loud- 
ness ratio judgments associated with simple 
numbers such as 2, 4%, 10, and the purposes of this 
investigation are adequately and appropriately 
met by direct tests so long as it can be shown 
that the results are consistent to a sufficiently 
high degree of statistical significance. This con- 
dition is found to be fulfilled by the results of the 
present investigation. It may be argued that even 
if the results are consistent they are not neces- 
sarily true. Of the various authors, GARNER [8] 
alone states that some of his observers could 
easily be persuaded to change their criteria and 
yet remain consistent. Careful observation and 
interrogation of the subjects who took part in the 
present experiments lend no support to this. The 
explanation of discrepancies between the various 
results of different experimenters is likely to be 
found in differences of procedure and the objec- 
tive conditions of test rather than in mere 
“suggestibility”’ or “faulty numerical sense”’ of 
the subjects. 

There is, in principle, a method of deducing the 
numerical loudness scale without actually making 
judgments associated with numbers. The method, 
which has been mentioned by GARNER [6], has 
some theoretical interest and the possibility of 
pursuing it had been considered in the present 
connection. The basis of the method is the esti- 
mation of “equal-loudness ratios” and “equal 
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loudness intervals’’, that is, the estimation of 
equal geometrical and arithmetical intervals on 
the scale of sensation. The difficulty about this 
procedure lies in presenting the comparison- 
stimuli in a satisfactory manner, since there are 
three of them. It may also be doubted whether 
the subtle distinction between ordering three 
stimuli by equal ratios and by bisection would be 
sufficiently intelligible to ordinary observers, so 
that the logical advantage of the method might 
well be outweighed by the experimental hazards. 
It would be interesting to know whether results 
would be obtained similar to those from direct 
ratio judgments but it should not be assumed 
that this is certain. It was decided that, for the 
stated purposes of the present investigation, this 
method is not appropriate, because the user of the 
sone scale does not have occasion to consider 
whether three noises stand in geometrical or arith- 
metical loudness relations but is concerned with 
more immediate comparisons. 


3. Choice of psychophysical method 


Two testing procedures are available for experi- 
ments of the present kind, namely, the “constant 
stimulus” and “average error’? methods!. Each 
of these basic methods may have a number of 
variants of a detailed kind. Opinions differ as to 
which method yields the more reliable data in 
particular experiments. In a paper devoted to 
this topic Kerrtioce [10] concluded that the 
constant-stimulus procedure was more reliable 
for experiments in loudness. STEVENS [11], on the 
other hand, advocates the opposite but with the 
reservation that practical or theoretical consid- 
erations may modify the choice. CHURCHER and 
Kine [12] have also discussed the relative merits 
of the two procedures. 

In the present case, several factors weigh in 
favour of using the constant-stimulus method. 
The question involved from the acoustical engi- 
neer’s point of view is whether a reduced noise is 
substantially or inappreciably less than the 
original, and he has usually to exercise his judg- 
ment in circumstances where he has no control 
of the stimulus. It seems most appropriate to 
make the laboratory judgments under similar 
conditions. 


1 These terms mean roughly the following. In C.S. pairs 
of stimuli are presented and the subject’s sole function is to 
judge whether the relation between them stands on one side 
or other of an assigned criterion; for instance, “Is A louder 
than B?”, or “Is A more than twice as loud as B?” In 
A.E, the subject adjusts one stimulus until it stands in the 
appropriate relation to the other. Thus, in the first 
example, B would be adjusted to sound as loud as A; in 
the second example, twice as loud. 
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Some of the other considerations appropriate 
to loudness ratio tests may also be mentioned. 
The manipulation of the switch or switch and 
volume control which are required in average- 
error procedure necessarily diverts attention from 
the judgment. In these circumstances some ob- 
servers may take a short-cut to an apparently 
stable final value which they do not repeat on 
subsequent occasions [5]. Another objection to 
the average-error procedure is that it provides 
the experimenter with little evidence regarding 
the quality of the judgments, whereas in the 
constant stimulus method, a spurious judgment, 
such as may happen through a slip of the tongue, 
is immediately apparent to the experimenter. 
Furthermore, the comparison stimuli may be 
presented at random intervals of intensity, thus 
forcing the observer to give independent consid- 
eration to each judgment, and encouraging him 
to be more discriminating in the search for barely 
detectable differences. There is, however, one 
drawback, namely, that the method undoubtedly 
occupies more time. A large number of the judg- 
ments made and recorded express the trivial fact 
that the stimuli were far removed from the desired 
value. This is the price to be paid for the out- 
standing advantage of randomisation. 

In the present investigation “don’t know” 
judgments were ruled out, but the comparison 
was repeated if necessary. “About right” judg- 
ments were discouraged but were finally admitted 
if the subject remained non-committal. This 


rarely happened. 


4, Experimental plan 


The experiments to be described covered the 
intensity range of practical importance. Some 
work at rather low intensities was included for 
completeness; the upper limit of about 110 phon 
was to some extent determined by instrumental 
limitations. 

Three quite distinct types of sound were em- 
ployed. These were a pure tone of 1000 c/s, a 
thermionically-generated white noise, the spec- 
trum of which was approximately constant per 
cycle band width up to a cut-off frequency of 
6000 c/s, and finally two reproductions of engine 
noises. Referring to Table I, the series I, I] and V 
is approximately parallel with III, IV and VI, 
and the purpose of using different types of sounds 
was to check the overall consistency of the ratio 
determinations by appeal to a logical proposition. 
Suppose A, and A, are two signals of 1000 ¢/s 
tone such that A, is judged to be twice as loud as 
A,. If, also, B, and B, are two signals of noise 
respectively judged to be equal in loudness to 
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A, and A,, then the existence of a numerical 
scale of loudness would imply that B, should be 
twice as loud as B,. It is interesting, and impor- 
tant for the practical utility of such a scale, to 
obtain direct experimental verification. Determi- 
nations of the equivalent loudness of the white 
noise were used to express the data of III, IV 
and VI in terms of equivalent 1000 c/s loudness, 
that is to say, in phons, and the converted data 
were then compared with the results of I, II 
and V. A similar conversion was used for the 
results with reproduced engine noises (VII). The 
truth of the proposition was amply demonstrated. 

Further satisfactory evidence of the reality of 
the numerical loudness scale was obtained by 
comparing the data for 2-fold and 10-fold ratios. 

The experimental programme was spread over 
a period of rather more than one year, the tests 
being completed in a non-systematic order in 
order to avoid any systematic bias due to habitu- 
ation of the team. 

The duration of a testing session was about 
15 minutes during which two separate tests were 
usually conducted. 


5. Details of experimental procedure 


The observer was seated at a distance of one 
metre from a direct-radiator loudspeaker in a 
sound-absorbent room, in which the ambient 
noise level was below the auditory threshold. 


Table I 


Plan of experiments 


Series Sound 


Type of judgment 


Doubling loudness 
Halving loudness 
Ill Doubling loudness 
IV Halving loudness 
Vv 10-fold magnification and 


1000 c/s tone 
1000 c/s tone 


white noise 


white noise 


reduction 1000 c/s tone 
VI 10-fold magnification and 
reduction white noise 
Vil Doubling loudness reproductions 
of engine noises 
VIII | Investigation of the effect of varying the interval 


of silence between the stimuli in loudness ratio 
judgments. 

IX | Investigation of the effect on the results of em- 
ploying monaural earphone listening. 

X | Investigation of the limiting behaviour of the ratio 
scale near the absolute threshold. 


In addition, 10 determinations were made of the equi- 
valent loudness of the white noise at various intensity 
levels covering the range used in series III, IV and VI, and 
similar determinations for the engine noises at appropriate 
ntensity levels. 
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The sound field of the loudspeaker was constant 
to within +0.5 dB over distances of 20 cm 
transverse to the axis, and obeyed the inverse 
distance law along the axis. The greatest harmo- 
nic distortion tolerated was just under 1 %, which 
occurred when the intensity level of the 1000 ¢/s 
tone was 110 dB above 0.0002 dyne/em? at the 
listening position; at lower intensity levels the 
distortion was much less than this. 

The sounds for comparison were gated by an 
electronic switch having a time constant of 0.15 s 
and low distortion during transients. A mechani- 
cal sequence generator controlled the switch, and 
this had normal mark and space periods of one 
second each. A single test group consisted of four 
marks and four spaces, 8s in all (Fig. 1). For 
experiments VIII only, the spaces were increased 
to 3s,so that one test group occupied 16 s. 


0 2 4 6 8s 
HS 
Fig. 1. Unit test sequence. 
(a = sequence generator, b = sound intensity, 
0: +°1, 4:5; + 2+ secondsi= fixed stimulus, 2° "3, 
6°*+7, +++ seconds = variable stimulus.) 


The time during which the maximum sound 
intensity persisted was chosen as short as possible 
to avoid unnecessary prolongation of the tests. 
The duration was adequate to ensure that full 
loudness sensation was reached [13]. 

The choice of the space period is of some im- 
portance. As will be seen, the results are usually 
repeated quite closely when this is increased from 
1 to 3s. Indeed it is quite possible that the re- 
sults would be repeated if a very prolonged 
silence separated the stimuli although the dis- 
persion would no doubt be relatively great. In 
practice a period of months may be relevant, as 
for instance in the subjective assessment of the 
efficacy of noise reduction projects. At the other 
end of the scale, however, too short an interval is 
objectionable in that adaptation to the first 
stimulus may influence the apparent loudness of 
the second. It is believed that the combination 
of 1 second mark and space is the most suitable 
experimental arrangement which at the same 
time gives results representative of useful prac- 
tice, 

The electrical circuit included separate attenu- 
ators controlling the “reference” and “test” 
channels which were switched for final amplifica- 
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tion by contacts operated by the sequence genera- 
tor. The attenuator controlling the “test” sti- 
mulus had a range of 60 dB in 1 dB steps. The 
range over which the comparison stimulus was 
varied depended on the type of test, and was 
typically 35 or 40 dB for 2-fold ratio judgments. 
After a few initial judgments the stimuli were 
mainly confined to a narrower interval of 15 or 
20 dB with occasional “‘easy”’ groups to prevent 
the observer as far as possible from memorising 
an absolute intensity. A>test consisted of 25 to 
40 groups, and a verbal judgment was recorded 
by the experimenter after each. The intensity of 
the comparison stimulus was varied in a non- 
systematic way but attempts to employ tables of 
random numbers were not very successful owing 
to the inefficiency which resulted from the fact 
that the “equality” datum varied considerably 
from one observer to another. In practice a 
qualified randomisation was used, in which sti- 
muli were offered in a random sequence around 
a rough placing of the “equality” point obtained 
fairly early in the test. 

In ratio judgments, comparisons were not 
offered for judgment which lay on the wrong side 
of objective equality since it was found that this 
procedure was liable to be interpreted as a fault 
condition in which the order of presentation had 
been accidentally reversed. In all ratio judgments 
the reference stimulus was presented first? and 
was maintained constant throughout a test. 

A team of 18 male and 7 female subjects was 
employed whose ages ranged from 18 to 52 and 
averaged 28. Some care was taken in the initial 
briefing of the team to avoid any mention of 
decibels, phons or sones. The instructions were 
cast in the simplest possible form, the observer 
being asked merely to state whether the loudness 
of the second sound heard was greater or less than 
twice (half, etc.) the loudness of the first. 


6. Experimental results 


The average results for each of the ratio tests 
are given first in Table II in terms of the decibel 
interval between reference and compared sounds. 

Of the two measures of average shown those 
in the column headed P(y) are the preferred 
values. They refer to the form of probit analysis 
described by GuiLrorp [14] by which the best- 
fitting normal ogive curves are deduced with a 
preferential weighting for central values (see 


appendix A). 


2 GARNER reports measurements in which the order was 
deliberately inverted. This procedure seems to be open to 
criticism, since the stimulus to be estimated has been 
withdrawn before the standard of comparison is disclosed. 
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The other measure of average shown is the 
arithmetic mean, obtained from estimations of 
individual decibel intervals from the test records, 
examples of which are shown in Table III. The 
values adopted were estimated from the test 
records to the nearest ¥% dB. 


7. Determination of equivalent loudness of the 
white noise and engine noises 


In order to compare the results obtained in 
ratio judgments with the pure tone and the noises 
it was necessary to obtain accurate equivalent 
loudness data for the latter. The testing proce- 
dure was very similar to that used for the ratio 
judgments, the main differences being that a 
slightly larger team (30) was employed, and that 
“equality” judgments were positively excluded. 
The results were found to depend on the order of 
presentation of the noise and tone signals, the 
difference being 0.98 dB in the sense that pre- 
ferential judgments of “louder” were made in 
case of the second sound heard. The “order” 
difference is highly significant by statistical test. 
A mean was taken in practice and the results for 
the white noise are shown in Fig. 2. 


120 


60 Ly Spectrum of white noise 


. 


Phons (Subjective determination by 30 subjects) 


Energy per unit 


—— = band-width 


20 40 60 80 100 
—— Sound pressure level of white ‘noise 
[dB rel.0.0002 dyne/cm?] 
Fig. 2. Equivalent loudness of white noise. 
3 Mean for 30 observers, with 98 % confidence 
+ limits; 
@ Data obtained by the method described in ap- 
pendix B. 


An interesting variant of the “constant- 
stimulus’? method was used to determine the 
equivalent-loudness of white noise at the highest 
intensity. This is described in appendix B. 


Table III. Specimen test records for loudness doubling 
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8. Analysis of results of ratio-judgments 


8.1. Table II shows the results of all the ratio- 
judgment tests, with the terminating levels of 
each interval expressed both in decibels and 
phons. When the intervals are expressed in phons 
in accordance with the experimental equivalent 
loudness determinations, the results for the diffe- 
rent types of sound are in close agreement. This 
may be seen especially clearly in Fig. 3a and 3b’. 
It is important to note, however, that if the inter- 
vals are expressed in decibels alone, this concord- 
ance is largely lost. 


Interval between sounds having 


20 40 60 80 100 120 
—— Equivalent loudness of reference sound [phon] 


Fig. 3. Results of 2-fold loudness ratio estimations ex- 
pressed in phons (average of 25 observers). 
(a) Doubling judgments [ 98 % confidence limits; 
(b) Halving judgments | 98 % confidence limits. 


© 1000 c/s tone (free field), 

A white noise (free field), 

@ 1000 c/s tone (earphone), 

mw 1000 c/s tone (free field) 3 second spaces, 
A white noise (free field) 3 second spaces, 

x reproductions of engine noises (free field). 


The consistency between the results using 
1000 c/s tone and the noises justifies single smooth 
curves to represent the data, as shown in Fig.3a 
and 3b respectively. However, it is not possible 
to amalgamate the fractional and multiple data 
without offending both. They differ in the same 
sense as found by previous investigators [5], [8], 
namely, that the interval for doubling at low 
absolute intensities exceeds that for halving and 
the contrary at high intensities. The difference, 
which appears from Fig.3 to be considerable in 
relation to the dispersion of the results, is certain- 
ly much less so in relation to the diversity of data 
by various authors. A tentative explanation of 


3’ The 98 % confidence bands shown on Fig. 3a and 3b 
extend + 2.49 x standard error in phons about the mean 
values, 2.49 being the value of SruDENT’s t for 24 degrees 
of freedom. They indicate that the true value lies within the 
bands shown with probability 0.98. 


Sol 


the effect is that it represents a predilection for 
moderate levels of the adjustable stimulus super- 
imposed on and mildly biassing the judgment of 
pure ratio. 

Not only do the various types of sound give 
results agreeing as to average but they also ex- 
hibit a high degree of consistency on the part 
of individual observers. This is illustrated in 
Table IVa, which shows the correlation between 
the various observers’ average results for each of 
the four sets of tests I-IV. 


Table IVa 


Correlation between individual results for white noise 
and pure tone 


Experimental series correlated Correlation coefficient 


I against III 0.64** 
II against IV 0;79%* 
I, II against III, IV O Mila 


** indicates that the significance level was 0.001 or less, 
that is, the chance of obtaining such values by pure 
chance is less than one in a thousand. 


Although, as mentioned above, the mean results 
for doubling and halving are significantly different, 
individual subjects nevertheless show substantial 
correlation in respect of this classification. This 


is illustrated in Table IVb. 


Table IVb 


Correlation between individual results for doubling 
and halving judgments 


Experimental series correlated Correlation coefficient 


I against IT 0.64** 
III against IV UO 
I, III against II, IV OxTe= 


** See note to Table [Va. 


Loudness scales for loudness reduction and 
magnification are shown in Fig. 4. These are 
derived straightforwardly from the smooth 
curves of Fig. 3. The curve standardised by the 
A.S.A. [2] is shown for comparison, the data 
having been arbitrarily aligned at 75 phon. It 
will be seen that there is substantial similarity 
between the shapes of the curves, especially if the 
A.S.A. curve is compared with the geometric 
mean of the other two. To bring the data into 
coincidence the logarithms of loudness must be 
transformed by a scale factor. 


8.2. Concordance of the results of 2-fold and 10-fold 
ratio judgments 

The data from experiments V and VI are illus- 

trated in Fig. 4. The general agreement between 
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the direct observation of 10-fold ratios and data 
read from Fig. 4 (which is based on 2-fold ratios) 
is very close. These comparisons are shown in 
Table V. The only significant discrepancy occurs 
in the case of 10-fold reduction from the loudest 
reference level (white noise), and this is believed 
to be due to 1-second periods between the stimuli 
being insufficient when the first sound is exces- 
sively loud. When the same experiment was 
repeated with 3-second intervals the discrepancy 
was brought below the level of statistical signifi- 
cance. The point is further discussed in section 8.3. 


200 
100 
50 


20 
10 


sone 


01 
105 
00. 0 


20 40 60 80 100 ~=—s-: 120 
—— phon 


Fig. 4. Loudness function relating phon and sone scales. 


this paper, loudness halving, 
—---- this paper, loudness doubling, 
x x American Standards Association data. 


(Note: The curves are arbitrarily aligned at 75 phon. 
The American Standard defines 1 sone as the 
loudness of a 1000 c/s tone, having loud- 
ness level 40 dB above threshold i.e. 40 phon 
for normal observers.) 


Results of 10-fold ratio tests 


@ 1000 c/s tone, 10-fold magnification (52.6), 

A white noise, 10-fold magnification (56.0), 

© 1000 ¢/s tone, 10-fold reduction (80.6) 

A white noise, 10-fold reduction (107.5 and 83.1). 


The figures in brackets are loudnesses of reference 
sounds in phon. The 10-fold reduction results shown 
thus [A] were obtained with the time sequence 
1 second mark—=3 seconds space. 


It should be noted that the agreement is be- 
tween 10-fold magnification and predictions from 
doubling: alternatively between 10-fold reduction 
and halving. The discrepancy between the halving 
and doubling intervals appears to persist in a 
general distinction between increasing and de- 
creasing sensations, perhaps for the reason sug- 
gested in section 8.1. 

The concordance illustrated by Table V_ is 
actually more deep-seated than the average 
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Table V 
Comparison of phon intervals for 10-fold loudness ratios 


Interval for 
10-fold loudness 
ratio in phons 


Type of sound and Error of 


equivalent loudness predic- 
of reference sound Predicted tion 
6 fi 
in phons doabkne Observed [phon] 


(Fig. 4) 


43.0 
41.3 


43.4 
39.6 


—0.4 
+1.7 


1000 c/s tone 52.6 
white noise 56 


Predicted 
from 
halving 
(Fig. 4) 


Observed 


1000 c/s tone 80.6 P +0.9 
white noise 83.1 110 +0.9 
white noise 107.5 + 8.8 
white noise* 83.1 times +3.1 
white noise* 107.5 | 1/10 + 3.3 


*In these cases only, the sequence of stimuli was 
l-second mark, 3-seconds space. 


figures show. By constructing curves of the type 
shown in Fig. 3 for each observer, individual pre- 
dictions of the 10-fold range were made from 
2-fold observational data. The correlation of the 
predicted and observed 10-fold intervals is shown 
in Table VI. It is strongly concluded that obser- 
vers are individually and collectively consistent 
in making ratio judgments over the wide sensa- 
tion range connoted by the numerals 2 and 10. 


Table VI 


Correlation of predicted and observed 10-fold ratio 
intervals 


Correlation coeffi- 
cient between pre- 
dicted and observ- 


Type of sound and 


; 3 Proba- 
reference intensity, 


in dB rel. ed dB interval for ee | 
0.0002 dyne/em? 10-fold loudness 
ratio 

1000 c/s tone 52.6 | times 0.003 
white noise 35 10 0.011 
1000 c/s tone 80.6 0.002 
white noise 60 0.001 
white noise 90 0.001 
white noise* 60 0.001 
white noise* 90 0.010 
All 10-fold tests 0.001 


* See note to Table V. 


8.3. The effect of increasing the silent interval 

It will be recalled from section 5 that some 
investigation was made of the extent to which 
the judgments of loudness ratio depended on the 
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time interval between the stimuli (experiment 
VIII). In all, five tests were made with a time 
sequence different from that used in most of the 
tests. The mark periods of 1 second were not 
changed, but the space was increased to 3 se- 
conds. These tests were suitably combined with 
others in which the normal 1:1 sequence was 
employed, in such a way as to provide controls on 
the effect under investigation. For ease of com- 
parison, the relevant data are shown in Table VII. 
Four of the five results, involving the ratios 2, % 
and 1/10 and different types of sound, show no 
significant difference for the altered time se- 
quence. The 2-fold results are seen to fit in quite 
well with the other data, as shown in Fig. 3. 

The two results involving the ratio 1/10 are 
shown in Fig. 4. Only one of these is in significant 
disagreement with the appropriate loudness curve 
based on %-fold data, namely, the case of 10-fold 
reduction from a disagreeably loud reference 
level. It seems possible that under this exceptional 
condition the shorter period of silence is not 
adequate for the ear to recover its sensitivity to 
the much quieter sound following. This is con- 
firmed by the much reduced disagreement which 
resulted when the same test was repeated with 
longer silent spaces. Probably because the second 
sound has not so marked a contrast with the very 
loud reference sound no such difference was ob- 
served between the results of halving judgments 
at comparable intensities when the time sequence 
was altered. 

It was concluded that with the one exception 
just discussed, the l-second space period was 
adequate. 


Table VII 


Comparison of results for normal and altered 
time sequence 


Type of sound and 
equivalent loudness 
of reference sound in 
phkons 


Interval in phons 


Type of test 


1s mark 
ls space 


1s mark 
3s space 


1000 c/stone 64.6 times % 12.6* 
white noise 62.4 times 2 14.5 
white noise 110.3 times % 12:37 
white noise 83.1 times 1/10 32.5 
white noise 107.5 times 1/10 30.9 


* interpolated from results of Series II. 
+ reference sound 111.9 phon. 


8.4. Classification of the results by sex and experi- 
ence 


In most of the ratio tests, the average results 
for males and females were about the same, but 
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the number of female subjects was too few to 
allow a precise conclusion to be drawn. 

As far as the acoustical experience of the ob- 
servers is concerned the conclusion reached is 
that this has no effect on the results. 17 out of the 
25 persons employed had no special knowledge 
of the subject or none at all, and had had no 
previous experience of this kind of test. 


8.5. The reliability of individual ratio judgments 
and the manner of variation of the 2-fold 
interval with absolute intensity level 


In addition to the concordance between 2- and 
10-fold ratio judgments which has already been 
discussed, there is substantial evidence that ob- 
servers maintain a constant criterion over a con- 
siderable period of time. 

An identical replication occurs in the case of 
the first two tests of series I, and in view of the 
null results of certain tests in series VIII and IX, 
it is admissible to regard these also as replications. 
The method of handling these data is outlined in 
appendix C and the conclusions may be briefly 
stated as follows: 


(i) the error standard deviation is of the order 
2.5-3.0 dB for 2-fold ratio tests, and greater 
though not proportionately greater for 10- 
fold tests, 


(ii) the results of individual observers show a 
strong constant tendency with regard to 
their departures from average, 


(iii) the trend of the curves shown in Fig. 3a and 
3b is broadly speaking the same for all 


observers. 


9. The effect of earphone listening 
(experiment IX) 


Comparison data for 1000 c/s loudness doubl- 
ing judgments at 40 dB using the free field and 
monaural earphone methods may be seen in 
Table II (cf. [IX and the third of I). In view of 
the null result only one such test was made. The 
result is also shown in Fig. 3a. The small differ- 
ence observed (0.85 phon) offers no basis for an 
explanation of the discrepant results in the 
literature. 

The earphone test was made using a pair of 
moving-coil telephones, only one of which was 
energised. The force of application was 400 g- 
weight. Probe microphone measurements showed 
that the sound pressure varied by +3 dB from 
one ear to another for constant excitation. 

Most observers commented unfavourably on 
the monaural earphone test, particularly because 
adventitious noises of physiological origin or from 
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the clothing seemed so much more noticeable 
than when listening to the loudspeaker. 


10. The nature of the loudness seale near the 
absolute threshold 


The limiting behaviour of the loudness scale 
close to the absolute threshold is of theoretical 
interest. Most of the previous data indicate a 
progressive diminution of the intensity interval 
for a ratio of half-loudness as intensity level de- 
creases, but the data appear to be inconclusive. 
If it were not for the essentially quantal nature 
of the sensation of loudness, the intensity inter- 
vals for a given ratio would ultimately have to 
contract to infinitesimal values to accommodate 
infinite subdivision of the scale. 


10 
Fig. 5. 
yA Diagrammatic form 
1 VA | of a possible loudness 
y function for an ob- 


server with normal 
threshold. (N. B. The 
width of the horizon- 
tal steps, represent- 
ing the differential 
intensity threshold 
and the particular 
form of the curve 
shown, are not in- 
tended to have any 
but diagrammatic si- 
gnificance.) 


sone 


40 60 
—— phon 


Fig. 5 shows diagrammatically a possible loud- 
ness function for an observer with normal audi- 
tory threshold. Suppose S has been found to 
sound twice as loud as P, and that Q is then used 
as the reference level for judgments of halving. 
Evidently both R and P will be “greater than 
half’’, whilst T is inaudible. In such a test, then, 
one would obtain a test record consisting exclu- 
sively of “greater” and “inaudible” judgments, 
no audible stimulus appearing to be “less than 
half”’. 

In the actual tests the 25 subjects were each 
asked to make loudness ratio judgments on 
stimuli-pairs, the first of which was adjusted to 
10 dB above the individual minimum audible 
field intensity at 1000 c/s (experiment X). In 
practice difficulty was experienced in preventing 
some of the subjects from making bisection judg- 
ments and the results are inconclusive. The mean 
phon interval for halving actually found was 
about 5, but this figure excludes 2 test records 
which did in fact consist exclusively of “‘ greater” 
and “inaudible” judgments. In view of the fact 
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that 5 phon is little more than the accepted 
value of the differential threshold [15], little 
reliance can be placed on the result. 


11. Assessment of previous data 


An attempt has been made to reassess the 
available data, which include three new determi- 
nations since the review given by CHURCHER [16] 
in 1935. Only such results are considered as are 
properly comparable, which entails the omission 
of a large part of the work of Larrp and of Ham. 
The data for comparison are shown in Fig. 6a-h 
together with an abridged annotation of the test 
conditions. 

In the case of Larrp, only the data for 1024 ¢/s 
are shown. Data for other frequencies are omitted. 
It may be noted that Larrp and his co-authors 
remarked on the consistency of their results for 
various frequencies, but this agreement would 
appear to be more apparent than real, since the 
results were compared in decibels rather than in 
terms of equivalent loudness. 

Ham and ParkInson’s procedure required the 
subject to estimate a percentage loudness relation 
between fixed stimuli. From the results given in 
their paper a loudness scale can be directly 
derived. The data for Fig. 6b were determined 
from such a scale, obtained by averaging their 
350, 1000 and 2500 c/s data. 

Fig. 6c and d shows the 2- and 4-fold data 
from GEIGER and FIRESTONE’s paper, the latter 
being included by appropriate adjustment. It 
is noteworthy that the 2- and 4-fold data are in 
reasonable agreement. It is also true that the 
10- and 100-fold data are concordant with the 
2- and 4-fold data as may be confirmed by 
constructing loudness scales from Fig. 6c and d 
and comparing the 10- and 100-fold results given 
in their paper. GEIGER and FIRESTONE seem to 
have over-emphasised the discrepancy between 
their results for different numerical ratios. There 
is some doubt, however, whether all the intervals 
reported by GEIGER and FIRESTONE are not too 
small. The author’s experience is that a 100-fold 
judgment could not be made from a starting-point 
as high as 55 dB, whilst a 1/100-fold judgment 
from 30 dB is virtually a determination of one 
difference limen or so above absolute threshold. 
The fact that GEIGER and FrrREsTONE criticised 
Lairp for presenting a rather wide range of 
stimuli suggests that they may themselves have 
tended too far in the opposite direction. 

The halving data from CuurcHER, Kine, and 
Davies’ paper are shown in Fig. 6e. Their 44-fold 
data are included by a suitable adjustment. 
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—— Phon interval for 2-fold ratio of loudness 


20 40 60 80 
— fquivalent loudness of reference sound 


Fig. 6. Previous determinations of phon interval for 2-fold loudness ratio. 


1700 phon 


The data for 1000 c/s doubling and halving 
given by Po.uAck are replotted directly. The 
data for noise have been converted to phons and 
phon intervals by means of Fig. 2 of PoLLAck’s 
paper, and are included in Fig. 6f and g. Although 
the data for noise and tone are somewhat dis- 
crepant a mean curve has been drawn to represent 
both, for halving and doubling respectively. 

The halving results of GARNER for a tone of 
1000 c/s are replotted in Fig. 6h. 

By averaging the ordinates of the curves in 
Fig. 6 the fairly smooth curves shown in Fig. 61 
result. There are six determinations of the halving 
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(a) Larrp, Taytor and WILLE 


Constant-stimulus method, 10 observers. Monaural ear- 
phone listening. Halving, 1024 c/s). 


(b) Ham and ParkINSON 


Method: Subjective assessment of fractional reduction 
as percentage. Free-field listening, ca. 50 observations 
per point selected from 175 observers. 


Data for 50 % reduction (derived from average loudness 
function for 350, 1000 and 2500 c/s). 


GEIGER and FIRESTONE 
Average-error method. 
Binaural earphone listening. 


ca. 40 observations per point (number of subjects not 
known); 1000 ¢/s. 


(e) CauRcHER, Kine and Davigs 


Average-error method. 34 observers. Binaural earphone 
listening. 800 c/s. 


@ Halving 
x 4-fold reduction (corrected) 


(f) Doubling 
@ 1000 c/s 


A white noise | PoLLack 


Average-error method. 7 observers. 


(g) Halving Monaural earphone listening. 


© 1000 e/s 
A white noise 


(h) GARNER 


Average-error method. 18 observers. Monaural ear- 
phone listening. Halving, 1000 c/s. 


Sequence: mark ls, space ls, mark Is, space 3s (repeat- 


ed). 


(i) Smoothed unweighted average of previous determi- 
nations. 


These data may be compared with Fig. 3. 


interval but only two of the doubling interval, 
contributing to these averages. The results of the 
present determinations conform quite well to the 
average halving curve but rather less well to the 
doubling curve (cf. Fig. 3). Unfortunately the 
only comparison data for loudness doubling are 
those from the papers by GEIGER and by POLLACK 
who happen to be those whose results for halving 
are much smaller than others. This may be the 
reason for the relatively poor agreement between 
the curves of Fig. 3b and 6i. 

The curves of Fig. 6i taken in conjunction 
with Fig. 3 suggest that a reasonable approxima- 
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tion to the function for loudness reduction is 
given by the formula logS=0.029(P—40) for 
the range 20-110 phon, the loudness being S in 
the sone scale and P in the phon scale. At lower 
intensities the function must be regarded as 
undetermined. 

It is difficult to find systematic procedural 
differences between the various results which cor- 
relate with the discrepancies, except possibly on 
the basis of the duration of the interval of silence 
between presentation of stimuli. For instance, 
Larrp, GARNER and the author, using | or 3 second 
intervals, obtain intervals higher than most, 
whereas Ham and Pouuack obtained low values. 
In Ham’s procedure there was no interval of 
silence; the same is inferred from POoLLAcK’s 
paper though the description given is not entirely 
clear on this point. In the view of the author, 
tests conducted without a silent interval are of 
doubtful relevance to the main purposes of the 
sone scale. 


12. Conclusions 


The conclusions drawn from this investigation 
of the relation between the sone and phon scales 
are as follows: 


(1) Highly consistent loudness scales are obtain- 
able by direct subjective ratio estimations. 
Estimates of the loudness scale obtained by 
loudness ratio judgments involving the nu- 
meral “2” agree closely with those obtained 
using “10”, 

(2) The scale for loudness reduction and also for 
magnification proved to be closely the same 
for entirely different types of sound, namely, 
a pure tone of 1000 c/s, a “white noise” with 
a spectrum level approximately constant up 
to 6000 c/s, and reproductions of engine 
noises. 


(3) There is a significant difference between the 
scales relating to increase and decrease of 
sensation respectively, especially at low levels 
of loudness. The interval on the phon scale 
for a given loudness ratio is approximately 
constant for decreasing sensation, but for loud- 
ness magnification it decreases as the absolute 
loudness is raised. 


(4) The shape of the function defining the loud- 
ness scale for individual observers is much the 
same for all observers. 


(5) The results of the present work combined 
with data from other determinations indicate 
that there is no appreciable difference be- 
tween the results for free-field and earphone 
listening. 
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(6) Judgments are not affected by a substantial 
change in the period of silence elapsing be- 
tween the test sounds if this is of the order of 
seconds. There is some indication from com- 
parisons with other work that the apparent 
loudness change is greater if no silent period 
is used. However, the use of judgments in 
which the sounds are presented without a gap 
does not seem to be relevant to the main ob- 
jectives of the sone scale. 


(7) The scale based on loudness reduction from 
the present investigation is in fair agreement 
with previous work. The results may be ap- 
proximated by the formula: 


log S= 0.029 (P—40) 
for the range 20-110 phons, 


P being the loudness level in phons and S the 
loudness sensation in sones. No simple for- 
mula expresses the results for loudness magni- 
fication. 
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APPENDIX A 
Estimates of the average in loudness ratio tests 


The method which was preferred was a form 
of probit analysis [14]. The data to be fitted were 
the total judgments cast in favour of the “test” 
sound exceeding the datum of comparison. In the 
present results, “about right” judgments were 
relatively infrequent and were counted as half a 
vote, thus avoiding the complication of deducing 
separate but barely distinguishable ogives for 
“exceeding” and “falling short” respectively. 
In every case the fit of the observed data to the 
best-fitting normal ogive ranged from good to 
exceptionally good by the chi-square test. Two 
specimens are shown in Fig. 7. The process em- 
ployed (URBAN’s constant process) estimates a 
median. The corresponding standard error, 
1.253s,///n, is given in Table IT. s,, which is an 
estimate of the standard deviation of the parent 
population, is inversely related to the abruptness 
of transition of the ogive curves. 
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Fig. 7. Specimen test results with ogive curves fitted by 
URBAN’s process. 


(a) White noise, 10-fold increase judgments, refer- 
ence intensity 35 dB rel. 0.0002 dyne/cm?. 


(Estimate of median intensity interval 
(50 % votes) = 41.3 dB.) 


(b) 1000 e/s tone, double-loudness judgments, re- 
ference intensity 69.6 dB rel. 0.0002 dyne/cm?. 


(Estimate of median intensity interval 
(50 % votes) = 14.6 dB.) 


The arithmetic mean result was obtained as 
described in the text and was distinct from, 
though generally quite close to, the median ob- 
tained from the ogive curves. The standard error 
of the mean, which has also been shown in Table IT 
is s/|/n, s, being the observed standard deviation 
of the individual estimates obtained from the test 
records. s, is generally more optimistic than the 
cognate statistic s,. 


APPENDIX B 
A variant of the “constant-stimulus” method 


A special method was applied to the determi- 
nation of the equivalent loudness of white noise 
at the highest intensity. Since this procedure has 
other applications it is considered worth de- 
scribing in detail. In the usual constant-stimulus 
method, it is necessary to have in hand quite a 
wide range, at least sufficient to carry the “test” 
signal well beyond the balance point for extreme 
observers. For loud levels of the reference signal 
this requirement becomes very difficult to meet. 
The case of 92.4 dB (white noise) affords a good 


illustration. Extrapolation from determinations 
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at lower intensities predicates a mean result of 
some 110 phon (Fig. 2) with a spread of perhaps 
20 phon. The extreme observer may be expected 
to judge equality at about 120 phon. Consequent- 
ly the circuit would have to be capable of repro- 
ducing 1000 c/s signals of up to 130 dB. This 
was beyond the undistorted capacity of the avail- 
able loudspeakers at the 1 metre listening posi- 
tion. 

A method is available for substantially reduc- 
ing the maximum level required; indeed, in prin- 
ciple, levels falling slightly short of the predicted 
average may suffice. In the present case that 
would be, say, 110 dB. Instead of presenting 
numerous pairs of stimuli to one observer, a 
single comparison repeated any number of times 
is offered. This is arranged to lie fairly close to 
the predicted average result. A single judgment 
is recorded for each subject. The process is 
repeated on subsequent occasions at one or two 
other intensities of the 1000 c/s signal in a fairly 
narrow range. The fractions of observers judging 
““test stimulus louder” are then used to deduce 
a normal ogive curve by the same process de- 
scribed in appendix A. The curve is interpolated 
(or even extrapolated) to the 50% fraction. This 
method has been found to be in excellent con- 
cordance with the conventional procedure (Fig. 2). 
The observed fractions and the best-fitting ogive 
are shown in Fig. 8 for the case cited. Some of the 
observers at no time judge the test stimulus 
louder, consequently the ogive curve must be 
assumed to be symmetrical. There is ample justi- 
fication for this by reference to the corresponding 
curves at lower intensities, which of course are 
determinate over their whole length. 


30 =100% re 
S § 20 —| 
=> 
Sc 
35 
2. 110.1 dB 
omg corresponds to 
s 2 10 50 % observers 
ee 
S" kao 
= 
ear Maximum sound | 
pressure level of test 
stimulus (112.6 dB) 


0 
100 110 120 


—= Sound pressure level of " test" stimulus (1000 c/s tone! 
[dB rel.0.0002 dyne/cm? J 
Fig. 8. Ogive curve determining the equivalent loudness of 


white noise. (Sound pressure level of white noise 
92.4 dB rel. 0.0002 dyne/cm?.) 
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APPENDIX C 
Statistical analysis of 2-fold test results 


The analyses of variance for series I-IV are 
shown in Table VIII. The components of variance 
between individuals and between reference inten- 
sities are significant in each case in relation to the 
error variance. Evidence of “observer < intens- 
ity’’ interaction is necessarily indirect. The resi- 
dual variances shown in Table VIII contain this 
interaction, if any, together with random error. 


Table VIII 


Analysis of variance and coefficients of concordance 


**Be- 
tween in- 

dividuals” 
variance 


**Be- 
tween in- 
tensities”’ 
variance 


Probabili- 
ficient |ty on the 
of con-| random 
hypo- 
thesis 


Residual 


variance 


(dB)? DF* | (4B)? DF | (dB)? DF 


I |} 9.5 24)96 6 | 7.9 144 | 0.593 | <0.0001 
I {5.5 24/42 6 | 6.6 144 | 0.530 | <0.0001 
III | 7.5 24])40 5 | 9.7 120 | 0.481 | <0.000i 
IV | 3.5 24/)3.9 4 | 4.3 96 | 0.538 | <0.0001 


* DF signifies the number of degrees of freedom. 


The further analysis of the residual variance 
requires an independent estimate to be made of 
the error variance. This is possible on the basis 
of sample replication data shown in Table IX. 


Table IX 
Replication data from doubling tests 


Estimates 
of error 
variance 


(dB)? 


Test replicated Known difference 


(i) 1000 e/s tone, 


x 2, 29.6 dB P3205 lapse of 3 months 
(ii) white noise, 
x 2, 40 dB 2.6 mark/space ratio 
(111) 1000 c/s tone, 
x 2, 39.6-40 dB 20.4 earphone vs free 
field plus lapse of 
12 months 
(iv) white noise, 
x %, 93-95 dB 10.9 mark/space ratio 


plus lapse of 5 
months 


The estimates of error variance in (i), (iii) and 
(iv) of Table [IX are homogeneous by BARLETT’s 
test. The estimate in (ii) is abnormally low, but 
can probably be left out of account in view of the 
much shorter lapse of time between the tests 
compared (a matter of minutes) as being too 
optimistic. The error variance is estimated to be 
of the order 14 dB? for 2-fold tests in general. 
This is actually rather larger than the residual 
variances shown in Table VIII, and must there- 
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fore account for the whole of it. The conclusion 
is that the estimates of error variance (Table IX) 
and the residual variances (Table VIII) are 
estimates of the same variance, and that con- 
sequently there is no significant component of 
variance arising from observer x intensity inter- 
action. 

The same conclusion is reached by calculating 
the coefficients of concordance [17] for the re- 
spective series of experiments. This statistic, 
which is obtained by placing the observed decibel 
intervals in rank order and summing over the 
various intensities, expresses the extent to which 
the observers giving high or low intensity inter- 
vals continue to do so at each absolute intensity. 
In the presence of a large observer = intensity 
interaction, the rankings would become confused 
and the coefficient of concordance would be low. 
The very high coefficients actually found for each 
experimental series I-IV strongly support the 
view that this interaction is not present to any 
extent, and that individuals on the whole main- 
tain a constant tendency with regard to their 
departure from average. The significance of the 
coefficients of concordance is estimated by the 
f-test against the variance to be expected with 
random ranking. The details are shown in 


Table VIII. 
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THE DIRECTIVITY OF SPHERICAL MICROPHONES 
by R. L. PRITCHARD 
Schenectady, N. Y., U.S.A. 


Summary 


By employing the reciprocity principle, alternative solutions for the response of spherical 
microphones considered recently by W. Kun are given in terms of the response for sound 
emission. An existing solution for the sound field radiated from a vibrating spherical cap is used 
to calculate the directional response of a cylindrical microphone mounted in a sphere. Similarly, 
the frequency response of a completely spherical (non-directional) microphone is calculated with 
the help of the electroacoustic reciprocity theorem. Results of numerical computations for each 
type of calculation are given and are compared with KuuHt’s results, which were obtained by 
direct integration of the sound pressure over the active surface of the sphere. 


Sommaire 


On indique, en appliquant le théoréme de réciprocité, deux solutions au probléme de la réponse 
a une émission sonore des microphones sphériques considéré récemment par W. Kuut. On calcule 
la réponse directive d’un microphone cylindrique monté dans une sphére, en utilisant des résultats 
antérieurs obtenus dans le calcul du champ sonore émis par une capsule sphérique qui vibre. De 
méme, on calcule au moyen du théoréme de réciprocité électro-acoustique la réponse en fréquence 
d’un microphone entiérement sphérique (non directif). On donne les résultats numériques obtenus 
pour chaque mode de calcul et on les compare avec les résultats de KUHL qui proviennent d’une 
intégration directe de la pression sonore sur la surface active de la sphére. 


Zusammenfassung 


Unter Verwendung des Reziprozitaétsprinzips werden weitere Losungen fiir die Richtcharak- 
teristik kugelférmiger Mikrophone, wie sie kiirzlich von W. Kunt behandelt wurden, aus der 
Richtcharakteristik bei Schallabstrahlung bestimmt. Eine Lésung fiir das Schallfeld bei Ab- 
strahlung von einer schwingenden Kugelkappe wird benutzt, um die Richtcharakteristik eines in 
eine Kugel eingebauten zylindrischen Mikrophons zu berechnen. In ahnlicher Weise wird mit 
Hilfe des elektroakustischen Reziprozitatstheorems die Frequenzkurve eines kugelférmigen 
Mikrophons ohne jegliche Richtwirkung bestimmt. Numerische Ergebnisse fiir jede Berechnungs- 
methode sind angegeben und werden mit den von Kuut durch direkte Integration itiber die 


wirksame Kugeloberflache gewonnenen Resultaten verglichen. 


In a recent paper [1], W. Kunt has calculated 
the directional response of a cylindrical micro- 
phone of diameter 2d mounted in a sphere of 
radius a as shown in Fig. 1. The method em- 
ployed was that of numerically integrating over 
the spherical cap 000, (which corresponds 
approximately to the circular area of the micro- 


Fig. 1, Cylindrical microphone of diameter 2d in a sphere 
of radius a; approximation by spherical cap of 
angle @,. 


phone diaphragm) the sound pressure existing on 
a rigid sphere [2] due to a plane wave incident 
upon the sphere from the direction @. In a similar 
manner, the frequency response of a completely 
spherical microphone also was calculated. 

It is the purpose of this paper to point out that 
exact solutions to the problems considered by 
Kuni may be obtained by employing the reci- 
procity principle. In the first case, the directional 
response of the spherical cap for reception of an 
incident plane wave is identical with the direc- 
tional response of a radiating spherical cap. For 
the second problem, the frequency response of the 
spherical microphone in sound reception may be 
calculated from the easily-determined frequency 
response for sound transmission with the help of 
the electroacoustic reciprocity theorem. 

A radiating spherical cap is defined as a rigid 
sphere having a radial velocity distribution that 
is a function only of the polar angle © and is 
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defined as 


v(a)= 1% V=C=23G; 
0 bs el SSS 


An exact solution for the sound field radiated 
from such a source with an arbitrary cap angle 0, 
has been given by Morse [3] and by STENzeEx [4]. 
The directional response f(9,0,)) may be defined 
as the ratio of the sound pressure at a large 
distance r in the direction 9 radiated by the 
spherical cap of angle O, to the pressure that 
would be produced at the distance r by a non- 
directional radiator of equivalent strength. (Equi- 
valent strength is defined as the product of the 
normal velocity vy and the radiating surface area 
of the cap.) In the notation of STENZEL, f(0,0,) 
may be written as 


1 
Mee) = T—c08 8) 0) 


= nt+1[Pn—1(c0s 95) —Pnii(cos Oo) 
Oe: U, (ka) +3V; (ka) | Pa(cos 6), 


n=0 


where k= 2z// is the phase constant (A= wave- 
length of radiated sound), P,(cos@) is the Le- 
gendre function of the first kind, P_,(cos@,) =1, 
and U,(ka) and V,(ka) are functions involving 
combinations of the Bessel functions Jm.y,(x) of 
half-integer order: 


Un (x) = (crx [2) 4 [x In ssie(x) —2 Ine (2) 


Va (x) =(—1)"*? (oex/2)" [e-Jn—n (+) 
+nJ_n—y, (x)]. 

By the reciprocity principle, it can be shown that 
f(9,0,) also represents the ratio of the average 
pressure existing over the surface of the spherical 
cap of angle O, due to a plane wave incident upon 
the sphere from the direction @ to the free-field sound 
pressure due to the incident plane wave. This latter 
ratio is the directional response considered by 
Kunu. A proof of the equivalence of the two 
directional responses is given in the Appendix. 

When @, is set equal to zero!, corresponding to 
a point radiator at the (north) pole of the sphere, 
the right-hand member of eq. (1) becomes iden- 
tical with the expression (eq. (35)) used by 
ScHwAnrz [2] for the relative pressure existing at 
a point on the sphere due to an incident plane 
wave. The reciprocal relation between the sound 
field due to a point radiator on a sphere and the 
sound pressure at a point on the sphere due to an 
incident plane (or spherical) wave was pointed 
out first by Ray.eicn [5], and has been noted by 
many writers [6]. 

"1 When OQ, is set equal to 0, the factor (1—cos@,)~! 
[P,, _, (cos Og) — PP, . , (cos Og)] reduces to (2n+1). 
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STENZEL has given many graphical illustrations 
of the variation of the function f (9, 0,) with @ for 
ka = 10 for different values of 0). However, the 


N 
3 


—— Directional response | f(@,20°| 


0 120° 


—+. @ 


180° 


Fig. 2. Magnitude of the directional response f(9,Q,) of a 
spherical cap of angle O, = 20° as a function of the 
angle © of arrival of an incident plane wave, for 
ka=10. Solid curve: tracing of Kunt’s caleulated 
curve; points: values calculated from eq.(1) above. 


particular case of 0)= 20°, which corresponds to 
the value considered by Kuut, is not included. 
On the other hand, numerical values of the co- 
efficients (the bracketed terms) in the series of 
eq. (1) have been given by STENzEL [7] for the 
case of ka = 10, 0, = 20°. Using these coefficients, 
the writer has calculated the directional response 
for this particular case. The results are indicated 
by points in Fig. 2, where the solid curve is a 
tracing of the curve calculated by Kuni and 
shown in his Fig. 6 as curve b. The agreement is 
quite good, except for values of O between 160° 
and 180°. In this region, the sound pressure at a 
point on the sphere changes quite rapidly with 0 
(cf. curve a of Kunt’s Fig. 6), and consequently 
reduced accuracy might be expected with the 
numerical integration used by Kuut. It might be 
noted that the values indicated by the points in 
Fig. 2 for O between 160° and 180° are in fairly 
good agreement with the experimentally-obtained 
values of Kuut (cf. the crosses, or curve ¢, in 
Kuut’s Fig. 6). 

In the conclusion of his paper, Kun. has ecal- 
culated the frequency response of a completely 
spherical microphone by graphical integration 
of the sound pressure over the surface of the 
sphere. This problem also may be solved exactly 
by employing the electroacoustic reciprocity 
theorem. The free-field voltage response (sensitiv- 
ity) of a microphone, or the ratio of the open- 
circuit voltage e.- at the microphone terminals to 
the free-field pressure p, at the location of the 
microphone, may be calculated in terms of the 
transmitting current response, or the ratio of the 
sound pressure p, produced at a large distance r 
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from the microphone acting as a sound source 
when a current 1, is applied to the terminals, by 
the relation [8] 


|eec/Pol =| Ps/ts|(2rA/g¢).- (3) 

Here, oc is the characteristic impedance of the 

medium, and 4 is the wavelength of the radiated 
and received sound. 

To determine the ratio p,/i., the radial velocity 

vy of the sphere (assumed to be constant over the 

entire spherical surface) may be introduced; thus 


Ps/ts = (ps/Vo) (Vo/ts)- (4) 
The first term, due to sound-radiation considera- 
tions only, may be determined readily from the 
solution for the sound pressure due to a zero- 
order vibrating sphere of radius a as [9| 


(ps/v) = (gea/r) (ka) [1+ (ka)?y”. (5) 
On the other hand, the ratio (vy/i;) of eq. (4) 
depends upon the nature of the electromechanical 
system. For an electrostatic transducer, the veloc- 
ity is directly proportional to the a—c charge, 
i;/jw, where @ is the circular frequency (w= kc), 
and is inversely proportional to the mechanical 
impedance [10]. With a stiffness-controlled me- 
chanical system, i.e., for frequencies below the 
fundamental-resonance frequency of the dia- 
phragm, as assumed by Kuut, the mechanical 
impedance is inversely proportional to frequency. 
Hence, in this frequency range, v, becomes 
directly proportional to applied current, inde- 
pendent of frequency. Under these conditions, the 
microphone sensitivity is directly proportional 
to the function (p./v,)) and to wavelength. Thus, 
combining eqs. (3), (4), and (5), 


leoc [po] = (4era?) jvg/is|-[1 + (ka)2}*. (6) 


. 

| -5 
=~ 
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= 
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4 6 
—~ ka = 2m a/r 


Fig. 3. Magnitude of the relative sensitivity of a completely 
spherical microphone (stiffness-controlled), of ra- 
dius a, in decibels relative to low-frequency value, 
as a function of the parameter ka=2za//. Solid 
curve: calculated from eq.(6) above; points: values 
caleulated by Kun. 
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A curve of the frequency-dependent term of 
eq. (6), i.e., the term within the brackets, ex- 


. pressed in decibels relative to the low-frequency 


value of unity, is shown in Fig. 3 as a function 
of the parameter ka= 2za//. The points in Fig.3 
represent values calculated by KuHL, which were 
obtained from his Fig. 8. In general the agree- 
ment is fairly good. 

In conclusion, the writer would like to em- 
phasize that the reciprocity theorem is a useful 
tool for solving problems of sound reception and 
transmission. Quite often, as in the cases des- 
cribed above, it is easier to solve the sound- 
radiation problem than to solve the corresponding 
problem in sound reception. 


Appendix 


A proof of the equivalence of the directional 
response of the spherical cap for sound reception 
with that for sound transmission may be con- 
structed by the use of the Green’s function, along 
the lines employed by Foipy and Priaxorr [11] 
for a general proof of the electroacoustic reci- 
procity theorem. For further details of this 
method, the reader is referred to their paper [12]. 

Points in space will be described by a radius 
vector R from the origin of an arbitrary coordi- 
nate system, while a point on the surface of the 
sphere will be denoted by r. The Green’s function 
for the spherical surface S is denoted by G(R, R’). 
Physically, G(R, R’) represents the pressure at 
the point R duetoa point source of sound of unit 
stength at point R’ in the presence of the sur- 
face S when the latter is perfectly rigid. From the 
symmetry of the Green’s function, G(R, R’) = 
G(R’, R), there follows the reciprocity theorem 
for the sound field. 

When the spherical cap is to be used as a sound 
radiator, the pressure p(Ro) at a point R, may be 
written as [11], [12] 


p (Ry) = (joo /4er) f [ G(r, Ry) vn (r) dS. (A1) 


The integral is taken over the surface S of the 
sphere, and vp, (r) is the normal velocity on the 
surface S (positive normal outward from the sur- 
face). An e) time dependence is assumed. Over 


— 


the spherical cap (0<O <0, of Fig. 1) vn(r) =v, 
while vn(r) =0( elsewhere on the sphere. On the 
other hand, the pressure at the point R, due to a 
point source of strength v,Scap located at the 


centre R, of the sphere may be written as 


Pol R,) =(jogvg Scap/4er)e 
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where S.ap denotes the surface area of the cap. 
The directional response fr for sound trans- 
mission may be defined as the ratio p(R,) [po( Ro): 
which is a function of the direction defined by the 
vector Raa. for a constant R= jp R: 
thus 


fr=[f {GG Rj) vn (tr) dS] /[v) Scape ?*/R]. (A3) 


5 


For sound reception by the spherical cap, a point 
source of strength Q is located at the point R,. 
In this case, the pressure pz(r) _existing at a 
point 7 on the sphere is given by 


Pg (rT) = (joeQ/4c) G(r, Ro). 

If the point R, is at sufficient distance from the 
sphere, the wavefront due to the point source 
will be substantially plane at the sphere. The 
average pressure over the spherical cap may be 
written as an integral over the entire sphere by 
introducing the function v, (r) /v) used above to 
define the location of the cap in the arbitrary 
coordinate system. Thus, 


Pe = (1/0 Seer) /{Ps@)rn(r)aS 


as 


= (1/t%9 Sap) (jog Q/4) | JG (F, Ro) v» (7) dS. 


Free-field pressure py at the centre R. of the 
sphere is given by the relation 


Po (Re) = (jmeQ/4ax) © *"/R. 
The ratio Pag/po(R.) for a fixed magnitude R is 
the directional response fi of the spherical cap 
for sound reception (microphone) as defined by 


Kun; thus (A6) 
far = [(1 v9 Scop) J J G (F, Ro) v (7) AS]/[e*"/ RI. 


(A5) 


Comparison 


eqs. (A3) and (A6) shows that the directional 


of the right-hand members of 
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response for sound reception is equal to that for 
sound transmission, i.e., fu=fr. 

It might be noted that formally, both the 
radiated pressure and the average pressure on the 
cap in sound reception are calculated by means 
of the same integral, cf.eqs.(Al) and (A4), 
which is difficult to evaluate in practice. The 
simplicity offered by the use of the reciprocity 
principle to calculate the microphone directional 
response in terms of the expression for radiated 
pressure resides in the fact that the latter may be 
evaluated by the alternative and simpler method 


of solving a boundary-value problem. 
(Received Ist February, 1953.) 
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Summary 


The absorption (x) of ultrasonic waves has been measured by a radiation pressure method at 
5, 10 and 15 Mc/s. The variation of «/f? with frequency (f) is discussed in relation to relaxation 
theory. In esters, generally, a pronounced decrease of «/f? with frequency occurs, thus justifying 
to some extent the relaxation theory. 


Sommaire 

On a mesuré par la méthode de la pression de radiation l’absorption des ultrasons aux fréquences 
5, 10 et 15 MHz. 

On étudie la variation de «/f? avec la fréquence au point de vue de la théorie de la relaxation. 
Il y a généralement dans les esters une diminution marquée de «/f? avec la fréquence, ce qui 
justifie dans une certaine mesure la théorie de la relaxation. 


Zusammenfassung 

Die Absorption von Ultraschallwellen wurde bei 5, 10 und 15 MHz mit der Strahlungsdruck- 
methode gemessen. Die Abhangigkeit der GroBe «/f? von der Frequenz wird im Zusammenhang 
mit der Relaxationstheorie diskutiert. Bei Estern findet man im allgemeinen einen merklichen 
Abfall von «/f? mit der Frequenz und kann damit bis zu einem gewissen Grad die Relaxations- 


theorie rechtfertigen. 


1. Introduction 


The classical formula for the absorption of 
sound in liquids is not found adequate to explain 
the anomalous nature of absorption in benzene, 
carbon disulphide and some other non-associated 
polyatomic liquids. In most of these liquids, the 
absorption coefflicient.of sound waves varies as 
the square of the frequency in accordance with 
the classical prediction assuming the absorption 
to be due to viscosity [1] and conduction of 
heat [2], but the amount of absorption in all these 
liquids is found to be much greater than that 
predicted by the classical theory. This led many 
investigators to find a formula which will explain 
this excess absorption. Out of these various 
theories the relaxation theory as given by HErz- 
FELD [3], KNESER [4] and BAvER [5] is now well 
established. According to this theory, there exists 
a relaxation of structural perturbation and of 
some internal degrees of freedom of the molecules 
of the liquid which will give rise to an additional 
absorption other than that due to viscosity and 
thermal conductivity. This additional absorption 
is given by KNneEsER [4] as 


po aet(Z) ie 


where [im = mC/ fm, fmis the relaxation frequency 


and c the velocity of sound. At low frequencies 
Xx. 2 [Am 
Teomenlakes 


can be simplified into 


cnalpt a) (2) 


This equation can be used to find the relaxation 
frequency from a plot of 1/x and 1/f?. The above 
equations lead to the conclusion that when 
f > fm, «/f? falls with the frequency finally ap- 
proaching the classical value when f > frm. 

In most of the non-associated polyatomic 
liquids, this relaxation frequency is far above the 
experimentally obtainable range which makes it 
impossible to put this theory to test. RAPUANO [6] 
made measurements of absorption in a number 
of liquids between 75 and 280 Mc/s, but in none 
of them did he find a decrease of «/f* with fre- 
quency. Even in carbon disulphide he did not 
observe any fall in the value of «/f? up to 
280 Mc/s in spite of the fact that the relaxation 
frequency in carbon disulphide is 52.5 Mc/s which 
is well within the range of his experiments. A 
number of experiments were made in water [7], [8] 
and in no case a drop of «/f* with frequency is 
noticed as indicated by theory. 


=A (say). Then KNEsER’s equation 
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The experiments performed in acetic acid [7], 
[8] up to 100 Mc/s and copper acetate and ethyl 
acetate [9] up to 50 Mc/s give some weight to 
the relaxation theory. In these liquids, the value 
of «/f? falls with the frequency as demanded by 
the theory but the value of «/f? never reaches the 
classical one. It is probable that the value of 
x/f? may attain the classical value at a still 
higher frequency than the one reached by the 
experiment. 

In view of the limited data available to test the 
relaxation theory, we felt the need to determine 
the variation of «/f? with frequency in a number 
of organic liquids not investigated so far and this 
explains the purpose of the present paper. 


2. Experiment 


The radiation pressure method is used to de- 
termine the absorption coefficient. Every liquid 
is doubly distilled and used immediately after 
distillation. The experiment is performed at 5, 10 
and 15 Me/s. The value of «/f? is calculated for 
each frequency for all the liquids measured and is 
tabulated in Table I. Fig. 1 gives the curves 
relating «/f? with f for the liquids mentioned in 
Table I. 


Table I 


Variation of «/f? with frequency 


Classical 
——<————<—_—_————} value of 
5 Me/s |10 Me/s/15 Me/s} «/f*x 1016 


No. Liquid 


1. | Cyclohexanone : 
2.| n-Butyl acetate 27.0 | 2.27 | 0.92 | 0.57 | 0.100 
3.| n-Butyl alcohol 27.5 | 1.10 | 0.72 | 0.80 | 0.840 
.4.| Iso propyl alcohol | 26.5 | 0.92 | 0.90 | 0.92 | 0.460 
5.| Pyridine 27.0 | 5.10 | 3.94 | 3.60 | 0.900 
6.| Toluene 27.0 | 1.29 | 1.24 | 1.23 | 0.078 
7.| Cyclohexanol 28.0 | 5.00 | 5.10 | — 4.890 
8.| Capryl alcohol 27.0 | 2.07 | 1.95 | 1.90 _ 
9.| Diethyl malonate | 28.0 | 3.68 | 0.61 | 0.52 | 0.210 
10. | m-Xylene 27.0 | 0.89 | 0.80 | 0.72 | 0.079 
11.| Ethyl acetate 25.0 | 1.90 | 1.07 | 0.78 | 0.080 


3. Discussion 


The results given in Table I show that the 
value of «/f? for non-associated polyatomic 
liquids like toluene and m-xylene, remains ap- 
proximately constant as also in alcohols and the 
viscous liquids like cyclohexanol. In esters, there 
is a pronounced decrease in the value of x/f? 
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when the frequency is increased from 5 to 15 Me/s, 
the value falling at 15 Mc/s to nearly one-third 
of its value at 5 Mc/s. Our results in ethyl acetate 
agree very well that of Beyer and Smiru [9]. 
The value of «/f? is tending towards the classical 
value as the frequency increases in esters. 


Mc/s 30 


Fig. 1. Variation of o/f? with frequency for liquids in- 
dicated in Table I. 


It is probable that the value of «/f* in esters 
may reach the classical value at still higher fre- 
quencies. Our results in esters support to some 
extent the relaxation theory. It is not possible to 
trace a fall in the value of «/f? in non-associated 
polyatomic liquids like toluene and m-xylene as 
their relaxation frequencies lie far above the 


range in which we have measured. 
(Received 20th October, 1952.) 
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